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1 Introduction 

The most frequently occurring final states in high energy particle collisions contain hadronic 
jets. Because of their large production cross sections, jet observables can be measured to 
a high statistical accuracy, and are therefore ideal for precision studies. Examples include: 
the measurement of the strong coupling, a s , from jet rates and event shapes in electron- 
positron annihilation; the determination of the gluon parton distribution function (and 
also a s ) in deep inelastic lepton-hadron scattering into two plus one jets; the measurement 
of parton distributions in hadron-hadron scattering from single jet inclusive production 
and vector boson plus jet production. Often the relevant observables are measured with 
experimental precision of a few per cent or better. Thus, theoretical predictions with 
the same level of accuracy are necessary to fully exploit the physics potential of these 
measurements. This usually requires the computation of next-to-next-to-leading order 
(NNLO) corrections in perturbative QCD. 

The straightforward calculation of jet cross sections in QCD perturbation theory is 
however hampered by the presence of infrared singularities in the intermediate stages of 
the calculation, which must be treated consistently before any numerical computation may 
be performed. At next-to-leading order (NLO) accuracy, using a subtraction scheme to 
handle infrared divergences is the approach of choice. Exploiting the fact that the kinematic 
singularities of QCD matrix elements are universal, one builds process and observable 
independent counterterms that simultaneously cancel both the kinematic singularities in 
real-emission phase space integrals and the explicit e-poles in one-loop virtual corrections 
(here the use of dimensional regularisation in d = 4 — 2e dimensions is implied). 

At NNLO accuracy, the calculation of fully differential cross sections is a challenging 
problem, and various extensions of the subtraction method at NNLO have been proposed, 
see e.g. refs. [1-12]. In very broad terms, when setting up any subtraction algorithm, two 
quite distinct difficulties must be addressed. First, one must define subtraction terms that 
properly regularise the real-emission phase space integrals and second, one must combine 
the integrated form of these counterterms with the virtual contributions, so as to cancel 
the infrared divergences of the loop matrix elements. In a rigorous mathematical sense, the 
cancellation of both the kinematic singularities in the real-emission pieces and the explicit 
e-poles in the virtual pieces must be local. On the one hand, this means that the subtraction 
terms and the real-emission contributions must tend to the same value in d dimensions, in 
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all kinematic limits where the latter diverge. On the other hand, the cancellation of explicit 
e-poles between the integrated subtraction terms and the virtual contributions must take 
place point-wise in phase space. This in particular implies that it is possible to write the 
integrated counterterms in such a form that they can be explicitly combined with virtual 
contributions, before phase space integration. From a practical point of view, the full 
locality of the subtraction scheme is also important to insure good numerical efficiency of 
the algorithm. Finally, the construction should be universal, i.e. independent of the process 
and observable being considered. This avoids the need to tediously adapt the algorithm to 
every specific problem. 

However, defining universal subtraction terms that are completely local in the real- 
emission phase space is rather delicate, and there is very little freedom to define these in 
such a way, that their integration over the unresolved phase spaces becomes convenient. 
One way to address these difficulties is to use counterterms that are not fully local, but 
whose complete analytic integration is tractable. For example, the antenna subtraction 
method [8-11] builds the subtraction terms from so-called antenna functions. These are 
simple enough to be integrated analytically, but they do not reproduce azimuthal correla- 
tions in gluon splitting, and the cancellation of e poles in the real-virtual contributions is 
also nonlocal. Because of this, actual numerical computations with the antenna scheme, 
such as the calculation of total rates [13-16] and event shapes [17-21] in electron-positron 
annihilation, require the use of an auxiliary phase space slicing. Another option is to de- 
velop dedicated subtraction schemes that are applicable only to some specific processes, 
such as the production of colourless final states, vector bosons [22, 23] or the Higgs boson 
[24], in hadron collisions. Then, one may even propose to dispense with the subtraction 
method altogether, and adopt a strategy such as sector decomposition (see [25] and ref- 
erences therein), where the full Laurent expansions (in e) of the real-emission and virtual 
pieces are computed directly, and their finite pieces combined. 

Nevertheless, despite the subtleties, it is possible to define completely local countert- 
erms for real radiation, by first carefully matching the various QCD factorisation formulae 
for unresolved emission [4, 26], and then extending the expressions obtained over the full 
phase space [1-3]. Recall that in the scheme of refs. [1-4], the NNLO correction to a generic 
m-jet cross section with no coloured particles in the initial state (work towards an extension 
to hadron- initiated processes is presented in ref. [27]), 



^.NNLO 



/ d(J^ 2 J m +2 + / da™ +l J m+ i+ / dcr^ V J m , (1.1) 

Jm+2 J m+1 Jm 



i m+2 J m+1 

is rewritten as a sum of finite integrals 



NNLO / j^NNLO , / j NNLO , / j^NNLO /-. \ 

a = dcr m+2 + / dcr m+1 + / da m , (1.2) 

Jm+2 J m+1 Jm 



where 

i NNLO _ J j„RR 

aa m+2 — } aa m+2 J m+2 ~ a<T m+2 J m 



i RR,A^ j j RR,A^2 t 

0-C7" m +2 Jm+l ~ aa m +2 J'rn 



to' ^ 

^T+\° = { [da^i + jT AaT+t] Jm+i - [6*Z* 1 + ( / d <W' 2 Al ) A 1 J ™ } e=0 > ( L4 ) 
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and 



da. 



NNLO 



d£J m+2 



da. 



RR,A 12 
m+2 



+ 



A RV > A 1 , 



do-. 



m+2 



(1.5) 



are each integrable in four dimensions by construction. In eqs. (1.1)-(1.5) J n (n = m + 2, 
m + 1 and m) denotes the jet function which defines the infrared safe observable being 
calculated. All approximate cross sections appearing in eqs. (1.3)-(1.5) above have been 
defined explicitly in refs. [1-3]. To finish the definition of the scheme, we must compute 
once and for all the one- and two-particle integrals appearing in eqs. (1.4) and (1.5). In 
previous publications, we have shown that it is possible to adapt and employ well-known 
techniques of loop integration, such as integration-by-parts identities and solving of differ- 
ential equations [28, 29], the method of Mellin-Barnes (MB) representations with harmonic 
summation [30-34] and also sector decomposition [25], to compute the integrals that arise, 
analytically and numerically. Indeed, all one-particle integrals, denoted formally by L 
above, have been evaluated with these methods [35-37]. The actual computation of inte- 
grated subtraction terms leads to a large number of rather elaborate phase space integrals, 
however these can to be computed once and for all. 

In this paper, we compute the integral of the iterated singly-unresolved subtraction 
term, / 2 do^,'^ 12 , over the P nase space of the unresolved partons. We find that the inte- 
grated approximate cross section can be written as the product of the Born cross section for 
the production of m partons, times a new insertion operator (in colour space), ■ We use 
the method of MB representations, as developed in this context in ref. [37], to compute the 
integrals appearing in the various building blocks of the insertion operator. In several cases 
we find multi-dimensional MB integrals that are very difficult to compute fully analytically, 
and hence complete analytic expressions cannot be obtained at present. Nevertheless, in 
these cases direct numerical integration of the appropriate MB representations provides a 
fast and reliable way to obtain final results. We stress that for phenomenological appli- 
cations, this is all that is required, hence, we make no severe effort to compute analytic 
expressions beyond those that are trivial to derive. As a numerical check, all integrals are 
evaluated using sector decomposition as well. Thus, each integral in this paper is obtained 
by two independent computations. 

Since the paper is quite long and rather technical, readers mainly interested in under- 
standing the general structure of the results or in some applications are advised to first 
read sections 2, 3.1 and 5. Section 2 sets the notation, and in section 3.1 we present the 
final expression for the integrated iterated singly-unresolved approximate cross section in 
eq. (3.13) and the new insertion operator in eq. (3.14). These two equations are the 
main results of this paper. In section 5 we discuss some examples, specialising the general 
formulae to the case of two- and three-jet production processes. The explicit definitions of 
the integrated counterterms are then presented in sections 3.2 and 4. The technical details 
of computing the integrated subtraction terms are given in appendices. 
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2 Notation 



2.1 Matrix elements 

We consider processes with coloured particles (partons) in the final state, while the initial- 
state particles are colourless (typically electron-positron annihilation into hadrons). Any 
number of additional non-coloured final-state particles are allowed, too, but they will be 
suppressed in the notation. Resolved partons in the final state are labeled with letters 
chosen form the middle of the alphabet, i, j, k, I, . . . , while letters chosen form the end of 
the alphabet, r, s, t, . . . , denote unresolved final-state partons. 

We adopt the colour- and spin-state notation of ref. [38] . In this notation the amplitude 
for a scattering process involving the final-state momenta {p}, \A4({p})), is an abstract 
vector in colour and spin space, and its normalisation is fixed such that the squared am- 
plitude summed over colours and spins is 

\M\ 2 = (M\\M). (2.1) 

This matrix element has the following expansion in the number of loops: 

\M) = \M {0) ) + . . . , (2.2) 

where denotes the tree- level contribution and the dots stand for higher- loop contri- 

butions, which are not used in this paper. 

Colour interactions at QCD vertices are represented by associating colour charges Tj 
with the emission of a gluon from each parton i. In the colour-state notation, each vector 
\M) is a colour-singlet state, and colour conservation implies 

(Y,T^\M) = 0, (2.3) 

3 

where the sum over j extends over all the final state partons of the state vector \M) (recall 
that we are considering processes where the initial state is colourless), and the equation is 
valid order by order in the loop expansion of eq. (2.2). 

Using the colour-state notation, we define the two-parton colour-correlated squared 
tree amplitudes as 

\M$ k) \ 2 = {MW\Ti-T k \M®) (2.4) 
and similarly the four-parton colour-correlated squared tree amplitudes, 

\ M( Sku,i/ = {MWm-TtoTj-TiHM®) • (2.5) 

We will also use the following product notation to indicate the insertion of colour charge 
operators between (.M^l and 

\M {0) \ 2 ®Ti-T k = (M {Q) \T v T k \M {<S) ) , 

(2.6) 

^^^{Ti-T^Tj-Ti} = {M^\{Ti-T k ,Tj-Ti}\M^) . 
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The colour-charge algebra for the product ^ n {T i) n {T k ) n = Ti-T k is: 

T i -T k = T k -T l if i^k; T\ = C h . (2.7) 

Here Cf i is the quadratic Casimir operator in the representation of particle i and we have 
C q = C F = T R (iV c 2 -l)/N c = (7V C 2 - l)/(2iV c ) in the fundamental and C g = C A = 2 T R N C = 
N c in the adjoint representation, i.e. we are using the customary normalisation Tr = 1/2. 

We also use squared colour-charges with multiple indices, such as Tf r = Ct ir and 
T\ rs = C f irg . In such cases the multiple index denotes a single parton with flavour obtained 
using the flavour summation rules: odd/even number of quarks plus any number of gluons 
gives a quark/gluon, or explicitly for the relevant cases at NNLO: 

q + g = q, q + q = g, 9 + 9 = 9, 

(2.8) 

q + g + g = q, q + q + q = q, g + q + q = g, g + g + g = g- 

2.2 Cross sections 

In this paper we shall need to use only cross sections of producing n partons at tree level 
with n = m, the Born cross section, and n = m + 2, the so-called doubly-real correction. 
We have 

d4°)(M) =M £ d <MM)^lM 0) (M)l 2 , (2.9) 

where A/ - includes all QCD-independent factors and d(p n ({p}) is the d-dimensional phase 
space for n outgoing particles with momenta {p} = {pi, . . . ,p n } and total incoming mo- 
mentum Q, 

i=l ^ 71 ' ^ i=l ' 

The symbol ^{n} denotes summation over the different subprocesses and S{ n y is the Bose 
symmetry factor for identical particles in the final state. Then the Born cross section and 
the doubly-real correction are simply 

daf n ({p}) = da8>({p}) and d^a(W) = d^ 2 (W) ■ ( 2 - n ) 

The final result will also contain the phase space factor due to the integral over the 
(d — 3)-dimensional solid angle, which is included in the definition of the running coupling 
in the MS renormalisation scheme, 1 

S - f d(d - 3)n - ^ (212) 

St - J (2^-3 " r(i - 6) • (2 - 12) 



In the MS renormalisation scheme as often employed in the literature, the definition of the running 
coupling includes the factor S t = (47r) E e _E7E . The two definitions lead to the same expressions in a 
computation at the NLO accuracy. At NNLO these lead to slightly different bookkeeping of the IR and UV 
poles at intermediate steps of the computation, but the physical cross section of infrared safe observables 
is the same. Our definition leads to somewhat simpler expressions at the NNLO level. 
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2.3 Momentum mappings and phase space factorisation 

The iterated subtraction terms are written in terms of momenta obtained via various 
combinations of the basic collinear and soft mappings of ref. [1]: 

{P} m+ 2^{P} { *1^{P}L S ' R) , (2-13) 

where both — ^> and — ^> may label either a collinear or soft mapping. (In general, both R 
and S are multiple indices.) As the above notation suggests, the final set of m momenta 
are denoted by tildes, while hats indicate the intermediate set of m + 1 momenta. In 
kinematic expressions where only the label of a momentum is displayed (we shall discuss 
several examples below), the tilde and/or hat is inherited by the label, and we write for 
instance i , ir and irs , where the latter two label single momenta. However, since these 
mappings affect only the momenta, but not the colour and flavour (apart from the flavour 
summation rules of eq. (2.8)), we shall omit the hat and tilde from flavour and colour 
indices. 

We also use labels such as (ir) to denote a single momentum that is simply the sum 
of two momenta, pu r \ = pi + p r . 

Importantly, both the collinear and soft momentum mappings lead to an exact factori- 
sation of phase space as follows: 

Mn+l ^ {P} { n r) ■■ # n+1 ({p} n+1 ; Q) = d(/) n {{p}^ , Q) [dp£J (p r , p ir ; Q)) , (2.14) 

Wn+i ^ {P} { : } ■ d^+xCMn+i; Q) = dM{P} { :\Q)[dp { ll( P r; Q)\ , (2.15) 

where here and below the bar denotes either a hat, when n = m + 1, or a tilde, when 
n = m. The one-particle factorized phase spaces can be written in the following form. For 
the collinear mapping we have 

[dp^{pr,Pir;Q)] = da ir (l-Q ir ) 2 ( ri - 1 ^ 1 - e )- 1 ^d0 2 ( K ,p r ;p (ir) )e(a ir )e(l-a ir ) , (2.16) 

where the pr ir \ appearing on the right hand side is understood to be written in terms of 
mapped momenta, that is pfcs = (1 — ai r )pf r + a^Q^ '. For the soft mapping we find 

[dpi r iCPr.;Q)] = dy rQ (l - y rQ )^ 1 ^~ 1 ^dMPr,K;Q)e(y rQ )e(l - y rQ ) , (2.17) 

where the momentum K is massive and K 2 = (1 — y r Q)Q' 2 - As the notation above indi- 
cates, ai r and y r Q will become integration variables, hence their precise definitions are not 
important and will not be recalled here. (See ref. [1] for details.) 

2.4 Kinematic variables 

Three types of kinematic variables are used to write the iterated subtraction terms. The 
precise definitions of these were given in ref. [1]. Here we recall only those formulae that 
are needed for defining every expression precisely in the present paper. 
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• Two-particle invariants, such as 

Si r = Zpi-Pr , S£ t? = 2p kt -p r , or s iQ = 2p v Q and s ik± = 2p i -k 1 _ . (2.18) 
The two-particle invariants scaled with Q 2 are denoted by yij = Sij/Q 2 . 

• Momentum fractions Zi >r and z^ ■? for the splittings p ir — > Pi+p r and p kt — > p k + Pti 

VlQ and z kt = - V%Q . » ( 2 - 19 ) 



ViQ+VrQ fe '* VkQ + y 



with z rt i = 1 — Zi tr and z^- j = 1 — zz t- Momentum fractions for three-particle 
splittings are denoted by 

Zfc.tr = ~^ 9 - l , (2-20) 

VkQ + VtQ + VrQ 

with the expressions for z r ^t and 2^ r fc obtained by cyclic permutation. In the fol- 
lowing, all momentum fractions will be integrated out, and so they will be expressed 
using the integration variables. 

We also use extensively the uncontracted and contracted eikonal factors: 

S^(r) = ^ , S jt (r) = g^(r) = . (2.21) 

Pj -p r p r -pi Sj r Sl r 

As mentioned above, the sum of two momenta is often abbreviated with the two indices 
in parenthesis, e.g. +Pr = P^ r y which is used systematically in other occurrences, such 
as 

• i c f,\ ^ s (ir)l Su + S r l /„ 00 ^ 

S(ir)k = Sik + Srk and Sr ir )i{t) = = ■ r — . (2.22) 

S(ir) t Slt {Su + SrtjSlt 

Finally, we express the integrated iterated subtraction terms as functions of the fol- 
lowing (combinations of) invariants: 

Xi = ViQ and Y^q = — — — . (2.23) 

UiQ UjQ 

In the centre-of-momentum frame (i.e. the rest frame of Q^), we find that Xi = 2Ei/y/s 
(s = Q 2 ) is simply the scaled energy of parton i, while 5y,Q = (1 — cosxij)/2, where Xij is 
the angle between momenta p^ and p^. 
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3 Integrating the iterated singly-unresolved approximate cross section 



3.1 The integrated approximate cross section and insertion operator 

We begin by recalling that the doubly-real emission cross section is defined precisely as in 
eq. (2.9), with n = m + 2. Then the iterated singly-unresolved approximate cross section 
times the jet function reads [1] 



d<w' 2 Al2 © Jm = m E ^ m+2 ({ P })-^— EfE^^iA^UM)! 



{m+2} 



s {m +2} V Lf^ 2 



^ k^t ' 



QJm{{p}), 
(3.1) 



where the notation J m means that the jet function multiplying the different terms in the 
sum depends on different sets of momenta. Explicitly, the three terms in eq. (3.1) are given 
by 



C^°U 2 |^ 2 (M)| 2 0J m ({p})^ 

C kt C^J m ({p}^ •»)) + C kt CS^J m ({p}^^) 



= E 

r^k,t 



c«c«rcs£? j m ({p}^.")) - c kt c rkt s^°\u{p} (£tM) ) 



+ E (^«cg2j m ({p}^.«))-c^^°v m ({p}^«)) 



(0,0) 



(3.2) 



5 t {O '°U 2 |^ 2 (M)| 2 0J m ({p}) 



,(0) 



(0,0) 



5 t C irt Gs£ 0) JmdpY 19 ^) ) - S t C trt S^J m ({p} {9 ' th 



5 t GS ir;t 5^ ) ^({p} (? ' t) ) + 5A^ r; ^^ ) J m ({p}( ? 



(0,0) 



+s t sir ) j m ({p}^)} 



(3.3) 



- 8 - 



and 



C k A m A 2 \MZ ! +2 ({P})\ Z ® J m ({p}) 



,(o) 



( fc r ,t)> 



+ E ( ^kAPSj? ^({P> (Y?,t) ) " C*^^^ J m ({p}< ? ''>) ) (3.4) 



(0,0) 



- C kt S t C krt S^J m ({p}^^) - C kt S t C rkt S^ 0) J m ({p}^^) 



(0,0) 



+ C kt S t S^ 0) J m ({p}^> t >) 



+ C fci S^°V m ({p}^). 



All momentum mappings in eqs. (3.2)-(3.4) lead to the factorisation of the original 
m + 2 parton phase space into an m parton phase space times two one-parton phase space 
measures, as discussed originally in ref. [1], and recalled in section 2.3 above. Symbolically, 
we may write 



d(j) m+2 ({p}) = d(j) m ({p})[dpi tm ][dp 



l,m+l\ 



(3.5) 



The jet function does not depend on the variables of the factorized one-parton measures, 
[dpi, m ][dpi im _|_i], so we can compute the integral of eq. (3.1) over the phase space of the two 
unresolved partons, independently of the jet function J m , that we shall omit in the follow- 
ing. The result of the integration is a long expression of kinematics-dependent functions 
- each corresponding to a specific iteration of unresolved limits of the squared matrix 
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elements — times colour factors, 



L 



da. 



' Al2 — TV £ d^ m ({p})-i 



m+2 



{m+2} 



{m+2} 



2k € Vq 2 



E | E 9 X] ( [ C fci C if!-]/ fc / t /r( T Lr) 2 + ^l C kt^kt;r}hl T kt T j T l 



~ [ C H C fcirC S i 2r]/fc/ i ( T 'L) 2 ~ [ C kt C rkt S kt ]frfkfti T r) 



+ E E \^V ht Cf^f k f t .j.f^ 
r^=k,t ij^r,k,t 

+ ^2[ c kt s u hk) t CAT o T i 

3,1 

+ EE f ^ ([StCi2] /i/r (T? r ) 2 + E[S.CSi°) ]^T 2 r T,T, 



-,2 \2 



3(0)1 „ 



+ E E [StSSV^^^.r^TyT,} + E l s t s rt] c ATiT k 

rj^t i,k,j,l i,k 



E E ( \^kt^t^frt]fkft( T krt) 2 ~ l^kSt^krt^H: ]( T k) 2 
k^t \-r^k,t ^ 

~ [ C A;A C rfct S fct K T r) 2 + y^[ C fct S t S rt ) ]^'^ r fc T 3 r i 



+ E E ( O [^fe*^t^iritl/i/» - T2 t T \ ~ [CfctStCS^-iS^jT 2 ^ j 

®IA<£?I 2 , 



r^k,t ij^r,k,t 



(3.6) 



where the operation (E) means insertion of the colour charges between (.Mrr?| an d 

see eq. (2.6). Three types of colour connections appear in eq. (3.6), and the functions on 

the right hand side — the "non flavour summed" (see below) integrated counterterms - 



take three different forms accordingly: 



(0,0) 



2 \ e 



X. 



(0,0) 



(2 \ £"i 2 
Q 2 



[xf\.TlTl\M^\\ 



ryWiOV)^! A/f (0) ,2 



A'. 



(0,0) 



(2 \ e 



(3.7) 
(3.8) 
(3.9) 



where e.g. [xj ^] represents a function that results in the integration of the counterterm 
X± . The quadratic Casimir operators that appear in eqs. (3.7) and (3.8) are factored out 



;o)]0'.0 



(together with [Xg ](^ fe )C?>0) dimensionless 



to make the integrals [X^ ]/■... and [X 
in colour-space. As the notation implies, ^ and [Xg may also carry flavour 

dependence. Incidentally, we note that for every integrated counterterm in eq. (3.6), we 
consider everything inside the square brackets to be simply part of the function's name. In 
particular, the lower indices inside square brackets loose their meaning. Nevertheless, we 
choose to keep these in order to exhibit from which counterterm each function derives. 

Eq. (3.6) is not yet in the form of an m-parton contribution times a factor. In order to 
obtain such a form, we must still perform summation over unresolved flavours, rewriting the 
symmetry factor of an m + 2-parton configuration to the symmetry factor of an m parton 
configuration. The complete details of this counting are not very difficult, but rather long, 
and are given in appendix A. As a result of these manipulations, we obtain functions - 
the flavour summed integrated counterterms — denoted by (X^ )) jf'^'", which are specific 
sums of the non flavour summed integrated subtraction terms, symbolically 



X (0) 



J/fe. 



(3.10) 



It is important to realise that the flavour indices on the left and right hand sides of the 
above equation need not match up. Indeed, the non flavour summed functions on the right 
hand side carry dependence on unresolved flavours, while the flavour summed functions 
on the left do not, by definition. Similar change in notation was introduced in the dipole 
subtraction scheme [38], where we find the definitions (see eqs. (7.21) and (7.22)) 



Vi(e) = V qg {e) 



= ^V gg (e)+n { V q q(e), 



if i = q,q, (3.11) 
if i = g, (3.12) 



where the functions Vj(e) on the left hand side represent flavour summed integrated coun- 
terterms, while the Vij(e) functions on the right hand side are not flavour summed. In the 
present paper, due to the extra complications of an NNLO subtraction scheme, we believe 
that it is helpful to make a sharper notational distinction between flavour summed and 
non flavour summed integrated counterterms. 
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After summation over unresolved flavours in eq. (3.6), we find that the final result can 
be written in the form 

I d<W2 Al2 = dal ® ig ({p} m ; e) , (3.13) 

where the insertion operator (in colour space) has three contributions according to the 
possible colour structures in eqs. (3.7)-(3.9): 



r(0) 

■ 12 



2vr 



(3.14) 



with /j denoting flavours, and C q = Cf = Tjj, C 9 = Ca = as in eq. (2.7). 
Eqs. (3.13) and (3.14) are the main results of this paper. 

In the following, we shall define and compute each term appearing in eq. (3.14). First, 
in terms of flavour summed integrated counterterms discussed above, we get: 



C 



(0) 

12,/, 



+ 



+ 



p(0) 



s 

cs 



(o),(j,0 

12 

(0),(j\I) 
12,/i 



5 (o),(t,fc)(j,0 

5 12 



p p(0) 
^kt^ktr 



q p(0) 



p p p;(0) 



q p /~q(0) 



p p q(0) 
^kt^rkt D kt 



fi 



c p q(0) 
°t^ J irt°rt 



Q P P5 c(") 
D t^irt KJ ~ > ir;t D rt 



p p(0) 
^kt^ir-kt 



fifk 



p q ret") 



p q(0)\ 



P pa(0) 



fifk 

(7.0 



p Q p(0) 



P P ps(0) 
^kt^ir-M^kt-^ 



r" ( CfcjSjCSj,, t S^ 



( Cfc t S t C fcri S^ 



fifk 
fifk 



fi 



(0) 



(i,0 



P c q(0) 
^kt^t^kt 



0\0 



P c q(0) 
^kt D t D rt 



a q (0) 



fi 

(7,0 
/< ' 

(*,fc)(j,0 



+ s*cs 



(o)\ 



(j.0 



Q P5 of") i 



(3,1) 
fi 



(3.15) 



(3.16) 
(3.17) 



(3.18) 



(3.19) 



On the right hand sides of these equations the flavour dependent functions are the flavour 
summed integrated counterterms discussed above. They depend on the kinematics through 
variables of the type a;, and ij^Q. The latter dependence derives from integrating an eikonal 
factor which is always multiplied with a colour-connected squared matrix element. In order 
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to make the results more transparent, we hid the arguments of the functions, but kept the 
relation to the colour-connected matrix elements, shown as upper indices. 

3.2 Flavour summed integrated counterterms 

Here we list the flavour summed integrated counterterms appearing on the right hand sides 
of eqs. (3.15)-(3.19), written in terms of the integrated subtraction terms. 

Collinear-type terms: 

1. Collinear-triple collinear: 



CktCktr) ~ i^kt^ktr\l99 + % l^kfiktrhgq + n i [^kt^ktrk'q'q > 

pk&ktr) = l^kt^At^ggg + n d^kt^itr\qq9 + 2n f l^kt^itr^aqq ■ 



2. Collinear-double collinear: 

M \ _ rp p(°) 



C (0) \ _ hn r (0) i 2r r r (0) i _ _ 

kt^irikt J — a l^kt^ir;kt\gg;gg + n f v^kt^irtktkqm 

' / qq 4 



+ (j^kt^ir-ktUqm + {C kt C^ kt ] ggm 



3. Collinear-soft-collinear: 

, (3.0 



c l «a<t) s tt ' , = i[c tt cs2t]g'» + n ( [c tl c S <°)itt') 



J kt^>kt;r ) „ ~~ 2 L fet kt;rigg "t L^fet^fctjrJ?? 



4. Collinear-triple collinear-soft-collinear: 



p p po(o) i rp n pc(o) i 

^kt^ktr^ktr ~ V^kt^ktr^kMriig > 



9 



5. Collinear-double collinear-soft-collinear: 

p pi 02(0) 1 rp p pq(o) 1 

^kt^ir-M^kt-s- ) , ~ V^kV^ir-.kt^kUrkg ' 

:(°) "\ _ 1 rp p 02(0)1 , _rp p ra(°) 



i.e. it is independent of the flavour /. 



(3.20) 



P pw \ _ rp pw i 

^kt^ir-kt 1 ~ y^kt^ir-ktmm > 

^fctCLfct) = ^[Cfc t C^ Jf^qg-gg + Tlf [C kt C - Aj , 

^fet^ir-fet) = 9 [Cfc^irLtlflS^S + n ^kt^\r;kt\iq\ig > (3.21) 



(3.22) 



(3.23) 



(3.24) 
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6. Collinear-double soft: 

, W) 1 



(c fc ,s£ ) ) WJ = ^C fet sS)]W) +nf[ c H sS!;)]^. (3.25) 



7. Collinear-triple collinear-double soft: 



^kfirkt^kt^j ~ 2^ kt ^ rkt ^kthgg + n f i^kt^rkt^kthqg ■ (3.26) 



Soft-type terms: 

1. Soft-triple collinear: 



q p(0A _ rc r {0h 



2. Soft-soft-collinear: 



(j,0 



0,0 1 
9 

3. Soft-triple collinear-soft-collinear: 



(s,<l) W) = ^[S t < ) J^ + n f [S,CB 



(0)1(3,0 
ir;tlqq 



9 

4. Soft-triple collinear-double soft: 



i.e. it is independent of the flavour /. 
7. Soft-double soft: 

(i,k)(j,l) 



fo o(0)\ W) _ [q cWiOV) 



(3.27) 



(3.28) 



(3.29) 



S tC irt S^) f -[S t C irt S%. (3.30) 



5. Soft-soft-collinear-double soft: 

^S t CSj r;t S^J = [S f CS iT .. t S^]^ , (3.31) 

i.e. it is independent of the flavour /. 

6. Soft-triple collinear-soft-collinear-double soft 

(SjC irt CS ir . ;i Srt^ = [S t Cj r . i CSj r;t S,|. t ' ) ] , (3.32) 



(3.33) 
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Soft-collinear-type terms: 



1. Soft-collinear-triple collinear: 



CfctSjC^J — [CjaStCjJ^gg + 2nf[C fct S t C Jfcr J M . 



(3.34) 



2. Soft-collinear-soft-collinear: 



n q ret") i \n q rc(o) i 



(CfctSjCS^^ — - [Ck t S t Cs!- r } t ]gg + n{[C kt S t (S\ r } t ] qq 



(3.35) 



fg 

i.e. it is independent of the flavour /. 
3. Soft-collinear-triple collinear-double soft: 



C fe AC fcr 4?) / - [C fct S t C fert S g>] + [C fct S t C rfct S £>] , (3.36) 



i.e. it is independent of the flavour /. 
4. Flavour-dependent soft-collinear-double soft: 

^°) > i w) = rc,.fi.fi( )iO-.0 



C«WJ = [C«W]W ; - (3.37) 



Actually, as shown here, and also seen in the precise definition of integrated flavour- 
dependent soft-collinear-double soft subtraction in eq. (4.37), the integral itself does 
not depend on the flavour. Distinguishing the flavour dependence serves book-keeping 
purposes: the flavour-dependent subtraction contributes to CS^M in eq. (3.18), while 

the flavour-independent one in eq. (3.39) contributes to S^'^'^ in eq. (3.17). 

5 . Soft-collinear-soft-collinear-double-soft : 

( CfctS t CS ir . t S rt ) = [C fet S t CSj r . t S^] , (3.38) 
i.e. it is independent of both flavours fi, fi- 

6. Flavour-independent soft-collinear-double soft: 

(c«S t sg>) (,,,) = [C«S t s£)]M. (3.39) 
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3.3 Analytic expressions to 0(e 2 ) 

In the next section we compute the functions on the right hand sides of eqs. (3.20)-(3.39) in 
terms of basic integrals that are calculated in appendices. Expanding these integrals in e we 
obtain the Laurent expansions of the functions on the left hand sides of eqs. (3.20)-(3.39). 
Analytic expressions for the expansion coefficients have been obtained to 0(e -2 ) accuracy 
in all cases, and we present these below. However, in the case of the single and double 
poles as well as the finite terms, we encountered several instances where obtaining complete 
analytic expressions was not feasible. This being the case, we made no severe effort to derive 
analytic expressions beyond those presented here. Higher order coefficients in the Laurent 
expansions will be given numerically in the form of a computer code elsewhere. In the 
following results we use d' = D' Q + d' x e (see appendix B and especially eq. (B.8)) and the 
abbreviations 

A,(n f ) llC A -4T R n f 



Bb(nf) 



C A 



N 



k=l 



3C A 

- {l-z) k 
k 



(3.40) 



(3.41) 



Collinear-type terms: 



1. Collinear-triple collinear: 

C A + 2C ft 



/ J I 



2C 



fi 



(3.42) 



2. Collinear-double collinear: 



' JiJk 



1 1 



2 (^ln Xi + In x k ^j - 3 S fit9 S fh , q - S (n f ) 8 fii9 5 fk , g 



- ( 3 + B (n f ) ) ( 8 fm 5 hig + S fi)9 6 fk>q 



+ 0(6- 2 ). 

(3.43) 



3. Collinear-soft-collinear: 



,(o) 



CfctCS^M (x u Y jltQ ) 

/ J i 



+ 



lnY ji,Q ~ \ ( 3( %,9 + So(nf) 



+ 0(e" 



(3.44) 



When e.g. i = j, the functions (c fet CSj^ r ) ( ) (x i: Y u> q) and (c kt CSf^ * ) (x^Y^q) are 

different (where i ^ j,l is understood) but up to this order, the functional dependence 
on the variables is the same. 
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4. Collinear-triple collinear-soft-collinear: 

(c kt C ktr ^ r ) f ( Xl ) = 1 -l(hx Xi + Z(y ,D' - 1)) 

+ ^l^ fi , q + B (n { )S fi , g "j+O(e 



(3.45) 



5. Collinear-double collinear-soft-collinear: 



(3.46) 



As implied by the notation, this function depends on both Xi and Y^n. However, the 
dependence on the latter vanishes up to this order. 



6. Collinear-double soft: 





1 1 



3 3 



? + ^ ^ In y j7 , Q + 4S(y , - 1) + 4 - ^("f ) 



+ 0(e" 2 ). 
(3.47) 



7. Collinear-triple collinear-double soft: 



Ck+C~kf S 



(o) 



' fi Xsfi 

This function is independent of the kinematics 
Soft-type terms: 
1. Soft-triple collinear: 



^-^Uv(y ,D> -l)-lB (n { )\ 



+ 0(e" 2 ) . (3.48) 



o r (o) 



fi 



C A + 2C 



fi 



2C 



f 



1 2 



1 1 



lnxi + S(y ,-Do' 



C A 



(3.49) 



+ 0(e" 2 ). 



2. Soft-soft-collinear: 



(StCSg)^ } (^,^ iQ ) = -1 + |(ln^ + Efa,,^ 



+ 



+ 0(e" 



(3.50) 



ln ^,Q " o 3^,9 + ^o("f) 



+ 0(e" 
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3. Soft-triple collinear-soft-collinear: 



1 2 



lnxi + E(y ,D' Q ) 
+ ^l( 3 ^ + ^o(n f )«5 /iig )+0( e 



(3.51) 



4. Soft-triple collinear-double soft: 



c p c(o) 

°t y - J irt D rt 



Ca + 2C fi 



+ 0(e" 2 ). (3.52) 



This function is independent of the kinematics. 
5. Soft-soft-collinear-double soft: 



c pa c(0) 



(j,0 

./', 



Oi^AQ) 



1 2 

^4 



+ ^lny i/iQ + 0(e^ 2 ) 

6. Soft-triple collinear-soft-collinear-double soft 

2 1 



{^t^irt^ir-^rt^j ( x i) ~ g 



£ 4 £ 3 



lnx i + 4E(y ,£)o) + 3S(yo,^-l) 



(3.53) 



+ 0( e ~ 2 ). 
(3.54) 



7. Soft-double soft: 



o o(0) 



(i,fc)(j,*) 



i (s(yo, £>{,) + £(y , #o - 1)) + 0(e~ 2 ) , 



c c(0) 
D t D rt 



(«,fc) 



+ i (z(y , D' ) + S(y 0) - 1)) + 0(e~ 2 ) . 



(3.55) 



The expansion for ^S t S^ 
eq. (D.54). 



is valid for the restricted kinematics discussed around 
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Soft-collinear-type terms: 

1. Soft-collinear-triple collinear: 



Ji 



Ca + C fi 



c 



fi 



1 2 



1 1 



lnxj + Y,(y ,D' Q ) 



C F 



C A 



(3.56) 



+ 0(e" 



2. Soft-collinear-soft-collinear: 

j(o) 



Cfc t S*CS \ r ' t ) (xfe) — -j — -3 



ln:r fc + £(y ,£>d) 

+ ^K 35/fc,9+ ^° (nfH/fc ' 9 ) +o(e ~ 2) 

3. Soft-collinear-triple collinear-double soft: 

(c fct S t C fcrt sS } )^ = I - l^dft,,^) + E(y ,D' - 1)) + 0( e ~ 2 ) 
This function is independent of the kinematics. 



(3.57) 



(3.58) 



4. Flavour-dependent soft-collinear-double soft: 



( 



Co o(0) 
kt a t a rt 



<Yfl,Q) 



+ 1 ( In F AQ + 2S(y ,^) + 2S(y ,^ " 1)) + 0(e' 



5. Soft-collinear-soft-collinear-double-soft: 



p 0. rc; q(") 

^kt^t^ir^rt 



1 2 



(3.59) 



(3.60) 



This function is independent of the kinematics. 
6. Flavour-independent soft-collinear-double soft: 



f r c c(0) 



(7.0 



(^,q) = ~ + ^ ( In fy,Q + 2S(y , ^) + 2S(y , ^ - 1)) + 0(e~ 2 ) . 

(3.61) 



Substituting the expansions in eqs. (3.42)-(3.61) into eqs. (3.15)-(3.19), we obtain the 
following explicit expressions for the kinematics dependent functions appearing on the right 
hand side of eq. (3.14). For C^y. we find 



+ 



C ft +C A 1 

2x 



S(y , D' - 1) - In Xi + \ U S fi , q + B 



>/«,<? 



(3.62) 



+ 0(6^). 
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The function f.f simply vanishes up to this order in the e-expansion 



C 12,f i f k ( x h x k> Y ik,Q) = 0(e z ) 



(3.63) 



The two-parton colour-correlated soft function, , is 



S l 3 m (Yji,Q) = ^( i lny AQ " E (yo,£>o) + £(</o,£>o - 1) + I ~ iMnt) 
+ 0(6" 2 ) . 

For CS^'ft we obtain (the Y ji,Q dependence vanishes up to this order) 



1 / 1 



3 3 



(3.64) 



CS 



(o),(j.O, 

12,/, 1 



41nx 4 -2S(yo,^o)-2S(y ,^o-l) 



3<5/ i , 9 +So(nf)*/ !1 , 



+ 0(6~ 2 ), 



(3.65) 



if i is distinct from both j and Z, while for e.g. i = j we have 



1 



7 3 

-lnx.j - 2£(y ,£ > o) - -£(yo,£>d - 1) 



^3^,, + B (nf)<5'/ i ,J 



+ 0(e" 



(3.66) 



Finally, the four-parton colour-correlated soft function, , reads 



j(0),(t,fc)(i,fc) 



; 12 



(**,q) = - ^ ( ^Yk,Q + s(yo, + s(i/ 0j ^ - 1) 

+ 0(6^), 



(3.67) 



if only two indices, say i and Zc are distinct, while 



,(0),(i,fc)(j,fc) 

'12 



1 / 1 



+ E(y , + S(y , - 1)) + 0(e" 2 ) , 



(3.68) 



for three distinct indices. Furthermore, the above expression is valid in the case of restricted 
kinematics of three hard partons. 



4 Integrated counterterms 

In this section we list the explicit definitions of the functions that appear on the right hand 
side of an equation among eqs. (3.20)-(3.39). 
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4.1 Integrated collinear-type counterterms 

1. Collinear-triple collinear: 



[c«c£] AA/p = (^Q^J^WtU 1 1 



Skt s kt'l 



1 p s.o.(0) / ^ , 2 \ (4.1) 

\ ktr I 

x f(a ,a k t,d(m : e))f(a ,aj: t9 ,d(m,e)) , 

where Pf °fljl are the strongly-ordered three-parton splitting kernels averaged over the 
spin states of the parent parton (see appendix D.l and especially eq. (D.27)). The sub- 
traction terms contain the spin-dependent kernels, that together with the corresponding 
kinematic variables can be found in ref. [1]. In appendix D.l we prove that the integrals 
of the spin-dependent kernels give the same contribution as those of the spin-averaged 
ones, therefore, we can use the latter when integrating the subtraction terms. 

The /(cko, a, d) functions, defined in eq. (B.6), represent simple modifications to the 
original subtraction scheme of ref. [1]. As discussed in appendix B in detail, these mod- 
ifications do not destroy the cancellation of singularities, but serve improved numerical 
control, efficiency and stability, and result in simpler, m independent, integrated coun- 
terterms. The rest of the counterterms are modified similarly by appropriate / functions, 
and below we shall simply include these factors without further comment. 

(i) 

The integrated counterterm is computed in appendix D.3. In terms of the functions X c 
(i = 1, 2 and 3) of eqs. (D.31)-(D.33) we find 

2 2 2 



r P r (0)i _ C fi 



kt 



c 



fk I 



, S^^,i C Si/. J • C /u/r,l I c' ) ( I B?; £ ' a »•*^'• , ) 



EE 

1=1 j=-\l=-l 



5 ht9 h f k f t h f ktfr ( x c 3) ( x lsir ' e > "o, d , 1)- 4 3) (x Wr ;e,a ,do, 2) 



(4.2) 



The various coefficients read: 



c (0) 


= 2, 


c (0) - 


-2, 


c (0) 


= 1-e, 


c (0) 


= o, 


c (0) 


= 0, 


c (0) - 

99,0 ~~ 


Tr 

Ca 1 


c (0) 


= b {0) 
u qq ' 


c (0) 

C 99,2 


= -6 (0) 
"99 


c (0) 


= 4, 


c (0) " 
c ggfi - 


-4, 


c (0) 
c ggA 




c (0) 

C 99,2 


= -&(°) 
u gg 



(4.4) 
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with c m,j = c 99,i' and finall y 

JP) _JP) _i d (o) _ d (o) _ ,(0) _ ,(0) _ ,(0) _ ,(0) _1 s 

2. Collinear-double collinear: 



rp p(0) , _ / 167T 2e 



2 \ 2 



2 kt - 1 - kt ir ir 

x /(«o, a kt , d(m, e))/(a , a^ f , d(m, e)) . 



(4.6) 

The integrated counterterm is computed in appendix D.4. In terms of the function X^ 
of eq. (D.39) we find 



2 2 

[CktC^ktlfkffjifr = Yl S tfkftjcfifrj 1 ^ ( x kt > X ^ ; e ' Q o> rf o; i, , (4.7) 

j=-lZ=-l 
with coefficients given in eq. (4.4). 

3. Collinear-soft-collinear: 

kw*«& = - (^« 2t ) 2 i[MSi[d« + J^j T (f (48) 

x f(a , a kt , d{m, e))f(y ,y ?Q , d'(m, e)) , 

where S~^(r ) is the eikonal factor defined in eq. (2.21) and pf^ are the spin-averaged 
two-parton Altarelli-Parisi splitting kernels (see eqs. (D.18)-(D.20)). 

(5) 

The integrated counterterm is computed in appendix D.5. In terms of the function Xq 
of eq. (D.53) we find 

2 

Vkt^M = E &,^ 5) (^^JT,Q.^H,Q.nH, (3 ; £ .«o, ( lo,i/o,4;i) 1 (4.9) 

i=-l 

with coefficients given in eq. (4.4). The integral X c with full kinematic dependence, as 
written above, first appears in computing NNLO corrections to four-jet production. 

4. Collinear-triple collinear-soft-collinear: 

V b e J J2 S ktr Z r.kt S ktT kt Ju 



3 fctf "r,kt ° kt - 1 kt 
x /(«o, a fct , d(m, e))f(y ,y?n, d'(m, e)) 



(4.10) 
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This integrated counterterm is also computed in appendix D.5. In terms of the function 
lf ] of eq. (D.70) we find 

2 

[C k t C ktr^kt;r\fkft = Yl C fkft,j X C \ x ktQ^ e ' Q 0, 4, 2/0, d' ] j) , (4.11) 

i=-i 

with coefficients given in eq. (4.4). 
5. Collinear-double collinear-soft-collinear: 

/[dpS][d^ +1 ]- 

ir *r'i ° kt ± kt 



\ J J2 b 7r Z r,7 bkt - L kt 

x f(a , a kt , d(m, e))f(y , y? Q , d'(m, e)) . 



(4.12) 

This integrated counterterm is also computed in appendix D.5. In terms of the function 



Xf ] of eq. (D.70) we find 



2 

[C w C iriJte CS2j p ] /fc/l = c$ ftjj l^\x TtQ ,Y irt Q -,e,a ,do,yo,d' -,j), (4.13) 

with coefficients given in eq. (4.4). 
6. Collinear-double soft: 

.2 \ 2 



MM = (^) /[dp£!][dpS +1 ] 



x ^(tt)^^(^} t (^ Ml ^,fc±, fclt ;e)|i/) (414) 
x /(«o, afci, e))f(y ,yg tQ ,d'(m, e)) , 
where S~ u ~(kt) is defined in eq. (2.21) and (mI-P/°/J 

f) is the spin-dependent Altarelli- 
Parisi splitting kernel for gluon splitting (see eq. (D.72)). 

(i) 

The integrated counterterm is computed in appendix D.6. In terms of the functions X c 



(i = 7, 8 and 9) of eqs. (D.84)-(D.86) we find 



i=-l 

(4.15) 

9 

+ h hh E ( X c «o, d ; 1) - 4 4) (e, q , d ; 2) 



with &^ * t and c^\ t & given in eqs. (4.3) and (4.4), respectively. 
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7. Collinear-triple collinear-double soft: 



16tt z 



2 c 



SktS 



P (0) (V , y, v ,h^> ^ML 



(0) Vk 



Tt \ za a T 



r \ kt,r * kt "ktQ , 

x /("o, afct, d(m, e))f(y , yg tQ , d'(m, e)) , 



(4.16) 



with the flavour-dependent constants 6^°^ given in eq. (4.3). Of course, only + 
ft = 9 gives a non- vanishing result (see eq. (3.25)). In obtaining the form (4.16) of 
the integrated counterterm, we exploited that the integrals of the two expressions in 
eq. (D.17) are equal. 

(s) 

The integrated counterterm is computed in appendix D.6. In terms of the functions X c 
and X^ of eqs. (D.85) and (D.86) we find 

" 2 

H 4!kr4 8) ( e >ao,do;i) 
-j'=-i 



[^kt^rkt^kt ]frfkft — TT - 

C .fr 



C A 



with coefficients given in eqs. (4.3) and (4.4) 
4.2 Integrated soft-type counterterms 

1. Soft-triple collinear: 



+ bfl (4 9) (e, ao, do; 1) - xf ^ «o, 4; 2) 



(4.17) 



[StCirtlfif; 



~s7 



Q ) I [ d Pi,m ] [ d ^l,m+l] rp2 PfifrgiZi^rt, Z Tj it, Zt.ir, Si r , Sit, S r t] e) 

x T 2 P fif r ( z f,i> e ^ ' d '( m ' e ))/( a o> u~ ? ,d(m, e)) , 

(4.18) 



(S) 



where the functions ij.j * are the soft limits of the triple-collinear splitting functions, 
introduced in ref. [4] (see eq. (E.2)), and we used T 2 irt = T\ r because the soft parton t 
can only be gluon, f t = g. 

This integrated counterterm is computed in appendix E.l. In terms of the functions 
X^ and Xf ] of eqs. (E.6) and (E.8) we find 



2 „(0) 



[s*ct°!] /l/r = y, 

k=-l 



fi fr , fe 
C fir 



(C/i + c f r ~ c hr) x s \ x ^^ a o, d o, i/o, do; A;) 



(4.19) 



+ 2C hr X s\ x Tr^ e ^ a o, d o,Vo, d' ; k) 
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with coefficients given in eq. (4.4). 
2. Soft-soft-collinear: 

16VT 2 



rc rc(°) i0'.0 



Q 2e ) /ICI'KUW) 



1 1 



j r it 



2 ± fifr^9,i- 



(4.20) 



x /(2/o, 2/tQ, d'(m, e))/(a , a ?f , d(m, e)) 



This integrated counterterm is computed in appendix E.l. In terms of the functions 
Z [ s ] and Z { s> of eqs. (E.ll) and (E.14) we find 



[S^S'tfo' = Yl c Mr,k [i 1 ~ d j(ir) ~ S l(ir))H ( x lf > Y jl,Q' e ' «o> d o, Vo, d' Q ] k) 

*=-l (4.21) 

with coefficients given in eq. (4.4). Note that j and I are always distinct. 
3. Soft-triple collinear-soft-collinear: 



[S t C irt CS[ r . t ]fif r 



x f{yo,ytQ,d'(m,e))f(a ,a~ 9 ,d(?n,e)) . 



(4.22) 



This integrated counterterm is computed in appendix E.l. In terms of the function 
of eq. (E.16) we find 



(5) 



[SA, t CSi% /r = £ c^ r>fc 4 B, (x-; ej a 0j do,Ito,4;fc), (4-23) 



.(0) T (5) 



fc=-l 



with coefficients given in eq. (4.4). 
4. Soft-triple collinear-double soft: 2 

16VT 2 



s f 



-Q 



2c 



^ « , r 



Ca (_ 



Sir 1 %r,it 



1 ££,rt i i Zijrt 



j r r ,i 



/(yo, ytQ,d'(m, e))f(y , y?Q, d'(m, e)) . 



(4.24) 



We note a harmless misprint in the definition of the subtraction term S t C irt S). t ' in eq. (7.38) of ref. [1], 
where in the last term of the square bracket ^4 was used as compared to here. Our definition of the 
momentum fractions gives the same for these ratios: — ii± — 
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The integrated counterterm is computed in appendix E.2. In terms of the functions Tg 
and 1 { J ] of eqs. (E.21) and (E.25) we find 

[S t C irt S% t = 2lf\e,y ,d' ) + ^±1^ ( e ,y ,d' ) . (4.25) 

fi 

5. Soft-soft-collinear-double soft: 

[s t cs ir ,s^ = - (^Vf / [dpS][d^ m+ i]^-^^iK*) 

x f(yo, VtQ, d'(m, e))f(yo,y?Q, d'(m, e)) . 

The integrated counterterm is computed in appendix E.2. In terms of the functions X s 
and Z< 9) of eqs. (E.35) and (E.40) we find 



[S,CS ir; 4?]W) = (l-5 j(ir) -S liir) )lP(Y llQ ;e,yo,d' ) 

+ (<*j(ir) +^i(ir))^s ) {xj,Y-rj Q ;e,y ,d' Q ) . 
Note that j and I are always distinct. 
6. Soft-triple collinear-soft-collinear-double soft: 

(0), _ / ^ n 2eY [ U Jrh U (t) , 2 Z 7t9 2 1 - g^r 



(4.27) 



S i r Z fJ S (ir)t z t,ir (4.28) 

x /(j/o, ytQ,d'(m, e))f(y , y ?Q , d'{m, e)) . 

This integrated counterterm is computed in appendix E.2. In terms of the function 1$°^ 
of eq. (E.44) we find 

[S t C irt (S ir]t S^] =4 10) Or 7 ;e,y ,d'o) ■ (4-29) 

7. Soft-double soft: 

= (^Q-) 2 ^[dpSS][dpSl +1 ]i5 ? ,(f )Sjl (t) 

x /(yo, 2/tQ, d'(m, e))f(y ,y ?Q , d'(m, e)) , (4.30) 
[S4?]^ } = - ^Q 2 ^ 2 ^[dpg][dp^ + i]^^(f )(5 ir (t) +S kr (t) -Sik(t)) 

x /(yo, 2/tQ, d'{m, e))f(y ,y ?Q , d'(m, e)) . (4.31) 

In this paper we do not discuss the case when i, j, k, and I are all distinct, which first 
appears in computing NNLO corrections to four-jet production. For the specific cases of 
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two and three hard partons in the final state, we compute the integrated counterterms 
in appendix E.2. In terms of the functions X^ 11 ^ and 2^ 12 ^ of eqs. (E.52) and (E.62) we 
find 

[sd?] mu ' k) = l4% k (yi~ k , Q ^~i3 >Q ^3~k, Q ^^X) , (4.32) 

and 

[S^?]^) = 2Z^%{Y~~ k Q - e, y , d' ) + 4;2^( y T^ Q ; ^ W>> 4) • (4-33) 

4.3 Integrated soft-collinear-type counterterms 

1. Soft-collinear-triple collinear: 

[C kt S t C krt ] fkft - {—Q j J[^ m ][d ^^ x] -——^^P Mr ^ £ ) 

x f{yo,ytQ,d'{m,e))f(a ,a k9 ,d(m,e)) , 

(4.34) 

where we used T krt = T kr because the soft parton t can only be gluon, f t = g. For the 
same reason T 2 kt /T 2 krt = T 2 k /T 2 kr , to be used in eq. (4.39). 

2. Soft-collinear-soft-collinear: 3 

\ / J2 s fct S^f 1 i r 

x f(yo,ytQ,d'(rn,e))f(a ,a^ ? ,d(rn,e)) . 



(4.35) 



3. Soft-collinear-triple collinear-double soft: 

2 \ 2 



x f(yo,ytQ, d'{m, e))f(y ,y ?Q , d'(m, e)) , 

,(0)1 _ ^ 1 6 7r2 (n 2eV fn (fe)irj (*) 1 2 2 Z M 



(4.36) 



x f(yo, ytQ, d'{m, e))f(y ,y kQ , d'(m, e)) . 

The integrals over [dp^] (or [dp 1 ^]) and [dpf m+1 ] factorize in the two equations above. 
Therefore, [C^C^?] = [C fet S t C 

rfet^i?] ( as seen the simple exchange of indices 
k f-> r) and the distinction of these functions is purely formal, and serves bookkeeping 
purposes only. 



3 We note a misprint in the definition of the subtraction term C kt S t CS^^ in eq. (7.46) of ref. [1]. The 
quadratic Casimir has to be changed from Tf to T\ (see eq. (7.22) of ref. [4]). 
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(4.37) 



4. Soft-collinear-double soft: 

[c kt s t s^ = - f^VY / [dpS][dp?Ui]^ ?T (f )AfM 

\ Je J J2 1 s kt %t,k 

x f(yo, ytQ,d'(m, e))/(y , y?Q, d'(m, e)) 

[c fct s4?] W) = - (^VY / [d^S][dp?Ui]kT(*)f IT 

x f(yo,ytQ,d'(m,e))f(y ,y^ Q ,d'(m,e)) , 

where the first equation above defines the flavour-dependent soft-collinear-double soft 
function, while the second equation gives the flavour-independent one. Again, the in- 
tegrals over [dpj*2] (or [dp^]) and [dp^| ri+1 ] factorize in both equations. Therefore, 
[C fet S^S^ = [C kt S t S^]^^ (as seen by the simple exchange of indices k r) and 
the distinction of these functions is formal, and serves bookkeeping purposes. 



5. Soft-collinear-soft-collinear-double soft: 

x f{yo, ytQ, d'(m, e))f(y ,y?Q, d'{m, e)) 



(4.38) 



The soft-collinear-type integrated counterterms are computed in appendix F. The cal- 
culation of the integrals in eqs. (4.34)-(4.38) is fairly straightforward because the integrals 
over the soft phase space measures decouple in all cases. Using the functions (i = 1, 
2, 3) computed in eqs. (F.2)-(F.4), we find: 



C 2 

[C w S t C^] A/t = c { ^ frJ l { ^{y rrQ ;e,ao,d ,yo,d' ;j) 

fkr j — ^ 



(4.39) 



[CfctS t CS^]/ 4 / r = c fif r! j T &\yk^Q^^ a o,do,yo,d' ;j), 
j'=-i 

[^kt^t^krt^rt ] = i ( ^kt^t < ^rkt^kt'\ = ^CS ( £ ' ^o) > 

[CfcASj?]^ = [C^sg^' 1 ) =2g ) (y ?7iQ ; 6, yo ,4), 

[CfcfSfCSj r;t S^] =1^(6,2/0,^0) • 

with coefficients given in eq. (4.4). 

This ends the definition of the integrals of the iterated subtraction terms, that can be 
used to construct the insertion operator as given by eqs. (3.14)-(3.38). The computa- 
tion of these integrals is presented in the appendices. 
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5 Insertion operator for two- and three-jet production 



In this section, we present illustrative numerical results for the insertion operator i| 2 of 
eq. (3.14), at various specific phase space points, for processes with at most three hard 
partons in the final state. 

While it is true that the most general collinear and/or soft configuration at NNLO 
accuracy involves four hard partons, if these are all in the final state, which is the subject 
of this paper, then the two-loop amplitudes needed in the doubly-virtual correction have 
at least six massless (or four massless and one massive) external legs. Results for such 
amplitudes are not foreseen in the near future, therefore, we restrict our discussion to 
computing NNLO corrections to two- and three-jet quantities. 

In an explicit computation of a jet cross section at NNLO, we require the expansion 
coefficients (up to and including the finite part) of the Laurent series in e of the insertion 
operator. In general, these expansion coefficients are functions of various kinematic vari- 
ables (and also parameters such as ao, yo, d,Q and d' Q ) which depend on the particular phase 
space point. One may either attempt to compute these functions analytically, or numeri- 
cally. The former is important as a matter of principle only. For practical purposes (from 
the point of phenomenology) the latter is sufficient. Indeed, the higher order expansion 
coefficients (starting form 0(e~ 2 )) of the results we will present were obtained numerically. 

5.1 Two-jet production 

Let us consider the process e + e~ — > 2 jets. The corresponding squared matrix element at 
tree level is \A4^\l q , 2^) | 2 , i.e. the quark carries label 1 and the antiquark label 2. Both 
the colour algebra and kinematics are trivial. Colour conservation implies 



and Cj x = Cf 2 = C-p. Hence, the insertion operator is a scalar in colour space. On the other 
hand, momentum conservation requires that the two final state momenta are back-to-back, 
i.e. in a properly oriented frame we have 



TiT 2 



C F 





(5.2) 



which implies 



Xi= X 2 = Y 12 ,Q = 1 . 



(5.3) 



The insertion operator eq. (3.14) becomes 



2tt e VQ 2 



) 



e-i 2 r 



) 



2C|(cf) 9 + C^-2CSgf 2 ) + 4S 



12 



(0),(1,2)(1,2) 



> 



(5.4) 



2CfCaS 



(0),(1,2) 



12 
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i 


-4 


-3 


-2 


-1 





r (o,i) 

X 12,2j 


6 


76 
3 


32.10 ±0.23 


-87.90 ± 0.66 


-554.5 ± 1.8 


T (x,i) 
X 12,2j 


-2 


27 
2 


-52.40 ±0.01 


-150.7 ±0.09 


-339.5 ± 0.43 


X 12,2j 





-1 


-6.332 ±0.001 


-17.66 ±0.008 


1.013 ±0.069 



Table 1. Coefficients of the Laurent expansion of the 2^i2,2j functions appearing in the insertion 
operator (pi , Vi j e ) i n the case of two-jet production. The numbers for i = —4, —3 are obtained by 
evaluating the appropriate analytic expressions. We used the parameters cto = yo = 1, do = d' = 3. 



with all arguments being equal to 1. Substituting the Laurent expansions of the kinematic 
functions, we obtain 

2 o 

(5.5) 



J r i2 > (Pi,P2;e) 
where [39] 



(2 \ e 
Q 2 



G F 2^ 6 l X 12,2j + XX 12,2j + V n f X 12,2j 



i=-4 



z = — , y 



(5.6) 

Of Of 

With this decomposition the Abelian case is obtained by setting Cp = I, x = 0, y = 1. We 
can compute the two leading terms in the e expansion analytically: 



4° \(pi,P2\ e ) 



2k e W 



e n 2 



1 



72 + 21x - 6yn f 
" 1 



(5.7) 



24(1 - x)E(yo, D' -l)- 12(2 - 3x)S(y , #c 



0J 



±0(e" 



The rest of the expansion coefficients are computed numerically. We present the results in 
table 1. To obtain these numbers, we used the specific values of «o = Ho = 1> do = d' = 3- 



Finally, we show the value of the complete insertion operator for the case of QCD with 
fif = 5 light flavours: 



83 



97.68 ±0.27 460.2 ±0.87 



1317. ±2.9) ±0(e 1 ). (5. 



I^( Pl , P2 ;e)- ^ i2e3 

In the above equation, the coefficients of 1/e 4 and 1/e 3 were computed by evaluating 
the appropriate analytic expressions. However, we have also computed these term with 
the same numerical algorithms that we used for computing the higher order expansion 
coefficients. It is then instructive to compare this numerical result to the analytic one. We 
find: 

r (o), v. 1-5 6.917 



1 12 ?(pi,P2;e)L 



+ 0(6-) 



e* e- 
1.498 ±0.0014 6.932 ±0.11 



±0(e- 2 ). 



(5.9) 
(5.10) 
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Comparing the exact and numerical results, we see first of all that the two match up to 
the uncertainty of the numerical computation, and second, that the error estimate on the 
next-to-leading pole is very conservative. 



5.2 Three-jet production 



Next, let us consider the process e + e~ — > 3 jets. The corresponding squared matrix element 
at tree level is \A4^ (lq, 2g, 3 5 )| 2 , i.e. the quark carries label 1, the antiquark label 2 and 
the gluon carries label 3. The colour algebra is again trivial. Colour conservation implies 



and C/j = Cf 2 = Cp, while C/ 3 = Ca- Thus, the insertion operator is again a scalar 
in colour space. On the other hand, the kinematics is no longer trivial, since the relative 
orientation of the three final state momenta are not fully fixed by momentum conservation. 
(Note that the insertion operator is independent of the overall event orientation with respect 
to the beam.) Since the three-particle phase space in d = 4 dimensions is 5-dimensional, 
but three of the independent variables just correspond to the three Euler angles needed to 
specify the overall orientation, we find that out of the six kinematic variables 



TiT 2 



Ca — 2Cf 
2 



T l Tz = T 2 Tz 



Ca 



2 




Xi = y iQ , and = Y ij)Q 



Vij 





OC 2 OC 'i 



only two are independent. Nevertheless, we choose not to fix the independent ones below, 
in order to better exhibit the structure of the insertion operator. 



The insertion operator eq. (3.14) becomes 



r (0) 

12 



/ 2 \ e ' 
Osn ( V_ 

2n e VQ 2 



[ Ci2, qg (xi,X3, Yu) + cf 2 ] gq (x 3 ,x h Y l3 ) 
+ cf 2 ] qg (x 2 ,x 3 ,Y 23 ) + cfl gq {x 3 ,x 2 ,Y 23 ) 



+ (c A - 2C F ) [c A (csf (^(xg,^) + s(°)'( 1 ' 2 )(y 12 ) 

+ c F (cs( )^ 2 )(x 1 ,y 12 ) + cs( ).( 1 ' 2 )(x 2 ,y 12 ) 

- c A [c A (csf^ 3 ) (x 3 , y 13 ) + s^ 1 ' 3 ) (y 13 
+ c F (csf ^ (x l5 y 3 ) + cs(°)-(M) (X2) y 13) 
+ c A (cs( )'( 2 ' 3 )(x3,y 2 3) + s(°)'( 2 ' 3 )(y 2 3 
+ c F (csf ^ (x!, y 23 ) + csf ^ (x 2 , y 23 ) 

+ 2C A (2C F - C A ) 



x S 



j(0),(l,2)(l,3) 



; 12 



(yi2,y 2 3,yi3) + si u 2 



(0),(1,2)(2,3) 



(yi2,yi3,y 



23 I 



,(0),(2,3)(1,2) 



(y 23) yi 3 ,yi2) + si u 2 



(0),(1,3)(1,2) 



(yi3,y 2 3,yi2) 



+ 2C|fsS2 ),(1 ' 3)(2,3) 



(y 3 , y 2 , y 23 ) + sg ),(2 ' 3)(1,3) (y 23 , y 2 , y 3 ; 



+ (8C 2 - 8C F C A + 2Ci)sS 2 } ' (1 ' 2)(1 ' 2) (y 12 ) 



2CifsS° 2 ) ' (1 ' 3)(1 ' 3) (yi3) + sg' (2 ' 3)(2 ' 3) (y 23 ; 



Substituting the Laurent expansions of the kinematic functions, we obtain 



(5.13) 



(2 \ e "t 2 
Q 2 



<3£ 



i=-4 



X I X 12,3j + xX 12,3j + x X 12,3j + 2/ n f X 12,3j + x 2/ n f X 12,3j 



(5.14) 
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We can compute the leading two terms in the e expansion analytically: 



1 12 ■ (Pl,P2,P3;e) 



2ir 



Q~ 



CU (<o + 2x + x 2 )\ + 



101 67 o 

12 H x H 

6 12 



13 3 

y yn f - -zj/nf - (4 - 4x)S(y , £>d - 1) 

(4 _ 6x _ x 2 )S(y , £>[,)- f 8 + x - ^ In y 12 



(5.15) 



4x+ -x 2 j(lny 13 +lny 2 3) ^ + 0(e 2 ) 

In order to obtain eq. (5.15), we used lnYy = In — In Xi — In . The rest of the expansion 
coefficients are computed numerically. For purposes of demonstration, below we present 
numerical results in three specific phase space points. Since the result is insensitive to 
overall orientation, we will always choose the event to lie in the x — y plane, with 
pointing in the positive y direction. In all cases we use the specific values of ao = yo = 1, 
do = d f = 3. 

Symmetric point. First, we consider the maximally symmetric configuration 



Pi 



u 
P2 



(5.16) 



Pa 



'I _JL_ _I o) 

1 _I 0) 

2^3 ' 6 ' ' ' 

which leads to the following values for the kinematic invariants = 1, 2, 3 and i ^ j): 

12 3 

Vij = 3 , Xi = 3 ' and Yi 'i = 4 ' ( 5-17 ) 

The coefficients of the Laurent expansion of the insertion operator in the symmetric phase 
space point are shown in table 2. 

Finally, we show the value of the complete insertion operator in the symmetric phase 
space point for the case of QCD with m = 5 light flavours: 



Ii2 (PijJfcjPa.e) = + + ' 

In this equation the coefficients of 1/e 4 and 1/e 3 were obtained by evaluating the appro- 
priate analytic expressions. As in the case of two-jet production, it is again instructive to 
compare the exact result with one obtained by numerical computation. We find: 



83 132.3 114.5 ±4.6 1142. ± 14. 



6150.±51.) +0(e 1 ). (5.18) 



i£Wft.;e)| 



i ( S(Pi, P 2,m;e)\ 



27.67 132.3 n . 2n 

— + — + )> 

27.658 ±0.027 132.4 ± l.£ 



+ 0(e- 



(5.19) 



(5.20) 



As before, the results match within the uncertainty of the numerical calculation, and the 
error estimate on the next-to-leading pole is again seen to be very conservative. We will 
reach similar conclusions in other examples as well. 
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i 


-4 


-3 


-2 


-1 





X 12,3j 


6 


34.12 


82.98 ± 0.25 


34.59 ± 0.71 


-543.8 ±2.2 


X 12,3j 


2 


9.721 


1.209 ± 0.52 


-142.2 ± 1.5 


-696.6 ± 4.8 


T (x 2 ,i) 
X 12,3j 


1 


6.497 


17.80 ±0.23 


15.87 ±0.79 


-47.93 ± 2.9 


X 12,3j 





13 
3 


-32.40 ± 0.007 


-127.9 ±0.03 


-355.2 ±0.20 


T (xy,i) 
x 12,3j 





3 
2 


-12.01 ±0.004 


-46.90 ± 0.02 


-104.1 ±0.16 



Table 2. Coefficients of the Laurent expansion of the %i2.3j functions appearing in the insertion 
operator (j>i,P2,P3i e) for three-jet production in the symmetric phase space point. The numbers 
for i = —4, —3 are obtained by evaluating the appropriate analytic expressions. We used the 
parameters ao = yo = 1, do = d' = 3. 

Collinear point. Next, we choose a configuration where (in the rest frame of Q^) we 
have a hierarchy of angles such that 

AP2,P-3.) « /.(pi,P2),£(pi,Ps) , (5.21) 
i.e. where momenta p^ and p^ are close to being collinear. Specifically, we set 

p v = ^ (0.49841 , 0. , 0.49841 , 0.) , 

v ^ = ^(0.120923 , 0.0240634 , -0.118505 , 0.), (5.22) 

p » = ^f s (0.380667 , -0.0240634 , -0.379905 , 0.) , 

which leads to the following values for the kinematic invariants (i, j = 1, 2, 3 and i ^ j): 

yi 2 = 0.238667 , y 13 = 0.758153 , y 23 = 0.003180 , 

x x = 0.99682 , x 2 = 0.241847 , x 3 = 0.761333 , (5.23) 

Y12 = 0.99 , Yi 3 = 0.999 , Y 23 = 0.0172697 . 

The coefficients of the Laurent expansion of the insertion operator in the collinear phase 
space point are shown in table 3. 

Finally, we show the value of the complete insertion operator in the collinear phase 
space point for the case of QCD with ?if = 5 light flavours: 

,(0), x 83 278.3 1601.3± 5.2 7084. ± 18. / m \ n , u . . 

I\ 2 \pi,P2,P3;e) = —+—+ (26690.±71.J+O(e 1 ). (5.24) 

We present next the comparison of the exact 1/e 4 and 1/e 3 pole coefficients given above 
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i 


-4 


-3 


-2 


-1 





X 12,3j 


6 


36.79 


106.0 ±0.23 


120.6 ±0.68 


-431.0 ±2.0 


X 12,3j 


2 


25.38 


143.6 ± 0.55 


537.3 ± 1.6 


1505. ± 5.3 


T (x 2 ,i) 
X 12,3j 


1 


15.24 


119.5 ±0.29 


660.5 ± 1.1 


2903. ± 4.9 


X 12,3j 





13 
3 


-31.30 ±0.007 


-121.7 ±0.03 


-346.0 ±0.18 


T (xy,i) 
x 12,3j 





3 
2 


-17.72 ±0.005 


-109.1 ±0.03 


-470.9 ±0.21 



Table 3. Coefficients of the Laurent expansion of the %i2.3j functions appearing in the insertion 
operator 1^ {pi,P2,Ps', e) for three-jet production in the collinear phase space point. The numbers 
for i = —4, —3 are obtained by evaluating the appropriate analytic expressions. We used the 
parameters ao = yo = 1, do = d' Q = 3. 



with ones computed numerically. We find: 



r (o) / \i 27.67 278.3 9 , „. 

I ( iJ(pi,P2,P3;e)\ A = + —3- + 0(e" 2 ) , (5.25) 

to), 27.658 ±0.027 278.4 ±2.1 n . 2 , . 
lfi(pi,P2,Pz; e)| N = ~ A + ~ 3 + 0(e~ 2 ) . (5.26) 

We note that the two results match up to the uncertainty of the numerical computation, 
and as in previous examples, the error estimate on the next-to-leading pole is seen to be 
very conservative. 

Soft point. Finally, we consider a configuration where (in the rest frame of Q 11 ) we have 
a hierarchy of energies such that 

E 3 ^E 1 ,E 2 , (5.27) 
i.e. where momentum is close to being soft. Specifically, we set 

p» = (0.480625 , 0. , 0.480625 , 0.), 

p% = yfe (0.487897 , -0.0308419 , -0.486921 , 0.), (5.28) 
p v = (0.0314778 , 0.0308419 , 0.00629557 , 0.), 
which leads to the following values for the kinematic invariants = 1, 2, 3 and i ^ j): 



yi2 


= 0.937044 , 


2/13 


= 0.024207, 


2/23 


= 0.038749 , 


x\ 


= 0.961251, 


x 2 


= 0.975794 , 


X3 


= 0.062956 , 


Y 12 


= 0.999, 


Y 13 


= 0.4, 


^23 


= 0.630768 . 



(5.29) 



The coefficients of the Laurent expansion of the insertion operator in the soft phase space 
point are shown in table 4. 
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i 


-4 


-3 


-2 


-1 





X 12,3j 


6 


25.85 


34.59 ±0.23 


-84.25 ±0.66 


-566.8 ± 1.9 


X 12,3j 


2 


27.79 


136.8 ± 0.52 


330.6 ± 1.4 


46.20 ± 4.5 


T (x 2 ,i) 
X 12,3j 


1 


21.02 


195.4 ±0.26 


1174. ±0.96 


5354. ±4.1 


X 12,3j 





13 
3 


-57.59 ± 0.009 


-405.2 ± 0.06 


-2119. ±0.34 


T (xy,i) 
x 12,3j 





3 
2 


-24.07 ± 0.006 


-194.7 ±0.05 


-1083. ±0.31 



Table 4. Coefficients of the Laurent expansion of the Ii2,sj functions appearing in the insertion 
operator 1^ {pi,P2,Pz] e) for three-jet production in the soft phase space point. The numbers for 
i = —4, —3 arc obtained by evaluating the appropriate analytic expressions. We used the parameters 
a = Vq = 1, d = d' = 3. 



Finally, we show the value of the complete insertion operator in the soft phase space 
point for the case of QCD with ri{ = 5 light flavours: 

r(0)/ . 83 320.6 1995. ±4.9 8928. ± 16. / \ n , u . 

J r i 2 ) (pi > P2,p 3 ; e) = ^4+^3- + ~2 " (32182.±61.J ±0(e 1 ) . (5.30) 

We finish by comparing the exact coefficients of 1/e 4 and 1/e 3 that appear above with the 
values obtained via numerical computation. We find: 

AO), m 27 - 667 320 - 56 -2s /r 

J r i2 ; (PiiP2,jJ 3 ;e)| A = — + — + °( e )> ( 5 - 31 ) 

,f ^ si 27.658 ±0.027 320.6 ± 1.9 nf 2 , . . 

^(pi.M.Pal e)| N = ^ + ^ + 0(e~ 2 ) . (5.32) 

Our conclusions are identical to those in the symmetric and collinear phase space point: the 
values match up to the numerical uncertainty and the error estimate on the next-to-leading 
pole is again shown to be very conservative. 

We finish by briefly commenting on the size of numerical uncertainties. The uncer- 
tainties relevant for phenomenology are those associated with the complete insertion 
operator, in various phase space points. However, the requirements in terms of precision 
are different for the pole coefficients and the finite part. 

On the one hand, the pole coefficients are only relevant for establishing the cancellation 
of all e-poles between the doubly-virtual cross section and various integrated subtraction 
terms. As stressed earlier, our subtraction scheme is fully local, hence this cancellation can 
be checked point by point in phase space for any specific process. From a practical point 
of view, it clearly suffices to demonstrate pole cancellation in a relatively small number of 
phase space points, thus the pole coefficients of have to be computed as precisely as 
feasible in a small set of points only. Because of this, the runtime of numerical integration is 
not an issue, and increased precision may be obtained simply by adjusting the parameters 
of the numerical integration to include more sampling points for each integral. 
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On the other hand, the precision requirement on the finite part of the insertion operator 
is essentially set by the relative uncertainty associated with the numerical phase space 
integration of the doubly-virtual contribution. This is not expected to be below the per 
mille level, hence from a practical point of view it is pointless to evaluate the finite term 
of the insertion operator with a precision much greater than this. In all cases discussed 
above, the relative uncertainty of the finite part of the insertion operator is at the per mille 
level already. 

6 Conclusions 

In this work, we have performed the integration of the iterated singly-unresolved approxi- 
mate cross section of the NNLO subtraction scheme of refs. [1—3]. The final result can be 
written as the product (in colour space) of the Born cross section times a newly defined 
insertion operator, I±2- The insertion operator depends on the colours, flavours and mo- 
menta of the final-state partons, and is an elaborate sum of many different terms, each 
corresponding to the integrated form of a specific iterated singly-unresolved subtraction 
term of ref. [1]. 

We have also explicitly evaluated all integrated subtraction terms which are necessary 
to assemble the insertion operator for processes involving at most three hard partons in 
the final state. The knowledge of these integrals (i.e. their Laurent expansions in e to 0(e) 
accuracy) is necessary in order to make the subtraction scheme an effective tool, and we 
have computed them once and for all. 

We have achieved this task by deriving Mellin-Barnes integral representations for all 
integrals under consideration. In principle, it is possible to evaluate all MB integrals via 
the residuum theorem, and in a subsequent step to obtain fully analytic expressions by 
performing the summation of nested sums over series of residua. However in practice, we 
have encountered several cases of higher order expansion coefficients, where the summation 
cannot be performed analytically with present methods. Therefore, in this paper, we have 
concentrated on the direct numerical evaluation of the MB integrals in the complex plane. 
All MB representations for both the numerical and, if available, the analytic expressions 
have been checked by an independent evaluation of the integrals using sector decomposition 
as in ref. [35]. We have found that all integrals contributing to the insertion operator are 
smooth functions of their variables (in the colloquial sense). For practical applications, 
this means that all integrals (in particular the finite in e contributions) can be given either 
in terms of interpolating tables or simple fitting functions, which can be computed once 
and for all. We leave this step for later work. Finally we want to stress again that 
the tables we have shown here are for demonstration purposes only, and obtaining high 
resolution interpolating tables needed for the computation of an actual cross section is 
straightforward. Increasing the accuracy of each entry is feasible as well, the best way of 
doing this is under investigation. 

The integrals discussed in this paper appear when integrating the subtraction terms 
that regularise the doubly-real NNLO correction to the jet cross section, see ref. [2]. The 
final step in finishing the definition of the subtraction scheme is the computation of the 
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integrated counterterms corresponding to the the doubly-unresolved approximate cross 
section (those labeled by A2 in ref. [1]). In that case, the analytic structure of the integrals 
is essentially the same as those studied in this paper, though a few are admittedly somewhat 
more cumbersome. Nevertheless, we are confident that the techniques of the present paper 
will also be applicable to the computation of these remaining contributions. 
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A Summation over unresolved flavours 

In this appendix, we discuss how to perform the summation over unobserved flavours 
in eq. (3.6). It turns out that we only need to consider three different cases explicitly, 
corresponding to sums involving two, three and four partons. All specific results are then 
easily obtained by appropriate substitutions. 

A.l Generic flavour sums 

Consider an m-parton configuration with ?nj quarks of flavour f,Tnj antiquarks of flavour 
/ and m g gluons. From this configuration we can obtain an (m + 2)-parton configuration 
in the following ways. 

1. Increasing the number of gluons by two, 




(A.l) 



2. Increasing the number of quarks and antiquarks of flavour / by one each, 




(A.2) 



3. Increasing the number of quarks and antiquarks by two each. The flavour of the two 
quarks may or may not be identical. We will refer to these two cases respectively as 



the 'equal flavour' (e.f.) and 'unequal flavour' (u.f.) configurations, 



2, rrif — > rrif + 2, 



m j 



rrif + 2 , 



m 



9 ~r m 9 2, 



SIS' 



m 



fir 



m 



SIS' 



m S/f + 1 . 



e.f. 
u.f. 



(A.3) 



where / and /' are understood to be different quark flavours. This case is relevant 
only for the doubly-collinear-type configuration, i.e. the sum involving four partons. 



The ratios of Bose symmetry factors for identical final state particles in the various cases 
are 



S. 



{m} 



D {rn+2} 



s 



(2) 

{m+2} 



1 



(m g + l)(m g + 2) ' 
1 

(rrif + l)(rrif + 1) ' 



(A.4) 
(A.5) 



and finally 



>{m} 



(3) 

{m+2} 



,5' 



m g (m g - 1) 



(rrif + 2)(m f + l)(m f - + 2)(m f - + 1) 



m g (m g - 1) 



(m/ + + + l)(mj, + 1) 



e.f. 



u.f. 



(A.6) 



Two-parton flavour sums. Consider a generic integrated counterterm pf^]f.'."\ which 
depends on two indices, k and t. This integrated counterterm may or may not depend 
on the corresponding parton flavours and it may or may not carry an upper index, as 
explained in detail below eq. (3.9). Examples would be e.g. [C^sj^ ]fff t an d [S t S^]^. In 
the latter case, we see a situation where no flavour index is displayed, since both r and t 
are constrained to be gluons. In what follows, we will discuss the most general case, when 
both flavour indices are explicit. 

Such terms necessarily appear in eq. (3.6) under a double sum: 



E 

{m+2} 



1 



s. 



{m+2} 



EM 



SkSt 



(A.7) 



In this configuration, we go from m to (m + 2) partons as in eqs. (A.l) and (A. 2). Then, 
since we are considering iterated singly-unresolved terms, both indices must correspond to 
unresolved partons, which means that they are either both gluons or a quark-antiquark pair. 
Now, decomposing the summation over t and k into sums in which the flavour (including 
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specific quark flavour) of each index is fixed, we find 

E EE^ki = 

{m+2} °{ m +2} t k+t 

= E s l : c (i { r } E z)[ x S ) ]ki <y /*.fl*/t.s (a. 

{m} ^ 5 {m+2} * fc# 



{m} 1™I ^ { , m+2} // i fc^J 



where Yip stands for the explicit summation over specific quark flavours. Performing the 
summation over t and k simply amounts to counting the number of ways in which we can 
assign the proper flavours to k and t in the appropriate (m + 2)-parton configuration: 

E E = #(/*)m+2 #(/*; k + t) m+2 (A.9) 
t k^t 

where #(/t) m +2 denotes the number of partons of flavour ft in the (m + 2)-parton config- 
uration, while k 7^ t) m+ 2 is the number of partons, different form t, of flavour in 
the (m + 2)-parton configuration. Note that t and are assumed to be distinguishable, 
which is the generic case. Clearly we have 

#{g)m+2 #(9; k / t) rn+2 = (m g + 2)(m g + 1) , 

(A.10) 

#(<?/)m+2 #(<?/; k ^ t) m+2 = (m f + l)(m / - + 1) . 

The case of ft = <jf and = is obtained by exploiting symmetry of this factor under 
permutations of indices. Then, using eqs. (A. 4) and (A. 5), we find 



E ^-EE^'&i = E ^-{[4?fc»+» t ([4?]^ + K 

m+2} *^{m+2} ( ^ {m} 3{m} I 



(O)i(..0 
kt iqq 



{m4 

(A.ll) 

In writing eq. (A.ll), we have used that [X^]^ and [-X^l^g do not depend on the specific 
quark flavour (as implied by the notation), and hence the summation Yip m e Q- (A. 8) may 
be performed, yielding the factor of rtf. 

Defining the flavour summed counterterm as 

E ^EE[4?]ki-E^rK , ) W > (A ' 12) 

{m+2} °{ m + 2 } t k# { m } °i m } 



we obtain 



3r(°)V'" ) - ry(°)l(-) -u n A _u ry(°)i(-) 



(A.13) 

where the second line follows, since [X^]^' 1 = [X^. 1 ]^'' in all cases we need to consider. 



ry(0)l(...) , o„ ry(0)l(-) 
(O)i(-) _ r y (0)i(...) 
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Three-parton flavour sums. Next, consider a generic integrated counterterm [X^,]^"-*, 
depending on three indices, k, t and r. Examples are e.g. [C kt C^ r ]f k f t f r and [S t CSj£. t ]jy?- 
As before, we will discuss the most general case, when all flavour indices are explicit. 
These terms always appear in eq. (3.6) under a triple sum: 



L ktrifkftfr ' 
{m+2} i m+z S t k^tr^k,t 



(A.14) 



In this configuration, we again go from m to (m + 2) partons as in eqs. (A.l) and (A. 2). 
Then, we decompose the summation over t, k and r into sums in which the flavour (in- 
cluding specific quark flavour) of each index is fixed. We obtain 



^2 c r nl ^2^2 ^2 [ X bl]f h f t fr 



S 

{m+2} i m + 2 } t k^tr^k,t 



E 

{m} 



1 S {m} x x - \ - [y (o),(...) 

5 {TO}5 (1) 2^,1^ 2^\- A ktr\f k f t f r 



'{m+2} t k^t r^k,t 



I 5. 



{m} s w s£U 



1 £££ £ [ X S]k/< 



Si- 



ft t k^t rjik,t 

6 fk,q f < 6 ft,q f ,$fr,g + S fk,q f > S ft,9 S Sr,q fl + 5fk,g 5 ft,q f ,$fr,q f , 

+ 6 fk,q fl 6 ft,q f ,$fr,g + $fk,q f >dft,g5f r ,q f/ + S fk,9 6 ft,q f ^fr,q f , 
■) 

Mr 



+ EEEEW 



/Y/ * k^tr^k,t 

6 fk,q f/ 6 ft,q f ,Sfr,qf + 6 fk,q fl 6 ft,q f ^fr,q fl + S Sk ! qf S St,q f ' S Sr,qf 

+ 6 fk,q f > 6 ft,q f ,dfr,qf + S fk,q f S St,qf S Sr,q f + S fk,qf S St,q fl S Sr,q f + (<?/ <"> 9/) 

(0)i(...) 
ktr\f k ftfr 



+EEEK l0,,t 



S Sk,qf S Suqf S Sr,qf + S Sk,qf S Suq f S Sr,qf + 6 fk,q f 5 ft,qf 6 fr,q f + (<?/ ^ If) 



(A.15) 



Next, we use the flavour summation rules to rewrite the summation over the unobserved 
indices k, t and r in the (m + 2)-parton configurations into a sum over a single index £;ir 
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in the m-parton configuration. We have 



— (A ' 16) 



t k^t r^k,t 



ktr 



where the notation is the same as in eq. (A. 9) and in particular #(fktr)m is the number of 
partons with flavour f^ tr in the m-parton configuration. Again, t, k and r are assumed to 
be distinguishable, which is the generic case. Then we have 



#(g/)m+2 #(g; k / t) m +2 #(gj T / k,t] 



m+2 



#(?/)m 

#(ff)-m+2 #(g; fc 7; #(gj r 7; fc, Qm+2 

#(<?)m 

#(g/')m+2 #(qf,k^ t) m+2 #{g;r^ k,t) m+2 

#(9)m 

#(g/0m+2 #(qf,k^ t) m+2 #{qf,r^ k,t) m+2 

#(?/)m 

#(g/)-m+2 fc ^ ^)m+2 r ^ fc, t) m+2 _ 

#(?/)m 



= (m s + 2)(m fl + l). 
(m g + 2)(m 9 + l). 
(m/' + l)(mj/ + 1) 



(A.17) 



(m// + l)(mp + 1) . 



(m/ + l){mj+ 1) . 



The rest of the cases are obtained by exploiting symmetry of this factor under permutations 
of indices. Then, using eqs. (A. 4) and (A. 5), we find 



fk ft fr 

_, U i rn.-\- /. V . , , . , , . 



{m+2} J { m + 2 > t kj=t r+k 



^ Ss m x ^ 
i m } 1 1 ktr 



ry(°)i(-) -i- rY-(°)i(-) -i- ry(°)i(-) o_ („ i^fyClit-) , ry(°)i(-) 

l^ktrlqgg ' l^ktrhqg ^ i^ktrlggq ~r V*f X M l^ktrlq'q'q ~r ^ktrlq'qq' 



Sqqq 



4. irC)i(-) 4_ rv(0)i(-) . ry(o)i(-) , rv(o)i(-) ^ , ry(°)i(-) 

L^fctr-lgg'g' I^fcir-W? L^fctrJfgg' ' r^-fctr-JggY j "•" l/Htr 

'(o)i(...) , rv(°)i(™)" 



6 fktr,q f + (?/ Qf) 



+ 



ry(°)l(-) i_ rv(0)l(-) , [y(0)l(-) , [y(0)l(-) 

l^ktrlggg ~r f \ L^fctrJooo ~r l^ktrlqgq "r L^fctrJflW " r V^ktAqqg 



■ fy(0)l(-) , fy(0)i(-A 
"•" l^ktriqgq KHtr Js<jg y 



fctnS 



(A.18) 



In obtaining eq. (A.18), we have used that whenever any of t, fc or r are (anti)quarks, 
[^fetrl/fc'/t/r does not depend on the specific quark flavour(s) (as implied by the notation), 
except that we have allowed for the possibility that the 'equal flavour' and 'unequal flavour' 
counterterms are different, e.g. [xj^ r ] q ^ / [^ktr\q'd'q^ wn i cn implies that the summations 
Y2f m e Q- (A. 15) may be performed, yielding the factors of (nf — 1) and rif. 
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Let us define the flavour summed counterterms as follows: 



£ - £ £ £ KSifcL = E st- £ ■ 

{m+2} » k^-t r=£k,t {m} 1 ' ^tr 



Then we find 



v(0)V"' ) _ [y(0)i(...) , fy(0)i(...) , r r (o)i(...) , frtr _ i^ryWiW , r 

^fctr J — L^fctrJggg l-^ktrlgqg "•" K^fctrJgw T V l f 1 M l/HirJg'g'g l^ktrlq'qq' 



(...) _ 

... Jggg 1 V"-ktr\gqg 1 ^''ktrlggq ' V'"I V ^K^fctrJg'g'g 1 ^'-ktrlq'qq' 

.(o)i(...) , ry (o),(...) . Iy (o),(...) , r^(o)i(»0 \ , rv(°)i(-) , rvC)i(-) 



i fyWlH , [yWll-J , fy^H-j , fyWll-J \ _i_ fyl u M-i _i_ ryWlv-v 
l^ktrlqq'q' "r L^fctrJg'g'g L^HtrJg'gg' "•" V^ktriqq' q' J ~T V^ktrkqq l^ktrlqqq 



, ry(0)l(-) 
l-^ ktr iqqq ' 

y(0)\ ( - } _ [y(0)l(...) , _ Av(0)i(...) , r^(0)i(...) , r y (0)i(...) , [y (0),(...) r y (0),(...) 

^fctr ^ — l-^ktrlggg "•" '*f ^^HirJggg ' L^fctrJgpg "•" l^ktrlgqq "r l^fctrJ?99 L^fctrJggg 



■ ry(°)i(-)\ 

y^ktrigqq J 



(A.20) 



Four-parton flavour sums. Consider finally a generic integrated counterterm [x£?L]f.'.''^ , 
which depends on four indices, fc, i, i and r. Two examples are [C^C^-J/j./^/^ and 

[C fet S t CS ir . t S^]. We will discuss the most general case, when all flavour indices are ex- 
plicit. 

These terms always appear in eq. (3.6) under a four-fold sum: 



E sr^EEE E [*&k/./,/,- < A - 21 > 



{m+2} { m+2 } t k^tr^k,ti^k, 



t.r 



In this configuration, we go from m to (m + 2) partons as in eqs. (A.1)-(A.3), i.e. case 3 
must also be considered. Decomposing the summation over t, k, i and r into sums in which 
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the flavour (including specific quark flavour) of each index is fixed, we find 



u \m+2\ 



ktirifkftfifr 



{m+2} 

= £ 

{m} 



t k^t r^k,ti^k,t,r 



1 s {m] (Q) , _j 

5, c (l) Z> 2-> Z> Z-> ^ktirifkftfifr 
'"'{m+2} * k^tr^k,ti^k,t,r 

S fk,qf S ft,9 S fi,q f ' S fr,9 + S fk,qf S fh9 S fi,9 S fr,q f > + S fk,9 S ft,q f S fi,qf S fr,9 

+ S fk,9 S ft,q f S fi,9 S fr,q fl + (if 9/) + 9/0 + (?/ *» 9/ > 9/' <-> 9/0 

S fk,qf S ft,9 S fi,9 S fr,9 + S fk,9 6 ft,q f 6 fi,g S fr,9 + S fk,9 S ft,9 S fi,q f S fr,9 
+ S fk,9 6 ft,9 S fi,9 S fr,q f + (9/ <-> 9/)) + $fk,9 S ft,g S fi,9 S f, ,9 



1 5. 



{m} 



{m} 



|m| '-'{m+2} /' * k^tr^k,ti^k,t,r 



(0) ,(...) 
ktir\f k ftfifr 



{ S fk,qf S ft,g S fi,q f > S fr,9f + S fk,qf S ft,9 S fi,q f > S fr,q f > + S fk,9 5 ft,qf S fi,q f > S fr,q f > 
+ S fk,9 S ft,qfSf t ,q f ,$fr,q f , + S fk,q f > S ft,q f > S fi,qf S fr,9 + S fk,qf S ft,q f > S fi,qf S fr,9 
+ S fk,q f > S ft,q f > 5 fi,9 S fr,q f + 5 f k ,q f/ S ft,q f ,$fi,g5fr,q f + (9/ 9/) 

+ S fk,9 S ft,9 S fi,q f > S fr,qf> + S fk,9 S ft,9 S fi,qf S fr,q f > + S fk,q f > 6 ft,q fl S fi,g S fr,g 

+ S fk,q f > S ft,q f/ S fi,9 S fr,9\ 



+ £ 



1 ^{m} VVW V V rr(°) i(-) 

Z^ Z^ Z^ Z^ Z^ Z^ ^ktirifkftfifr 



(A.22) 



{m} {m} ^{m+2} / /' * kfrr^k,t&k,t,r 

S fk,qf S ft,qf S fi,q f > S fr,q f/ + S fk,qf S ft,qf S fi,q f > S fr,q f/ + 6 fk,q f 6 ft,qf S f„q f > S fr,q fl 

+ S fk,qf S ft,qf S fi,q f ' S fr,q f/ 

Four-index subtraction terms only arise in conjunction with the double collinear limit and 
are always completely independent of the specific quark flavours. We have used these facts 
to write eq. (A.22) in the above form. First, since the pairs of indices k, t and i, r will 
always correspond to true singly-unresolved collinear limits, we have discarded all terms 
where this is not the case. In effect, we have dropped all terms where both k and t or both 
i and r are (anti)quarks. Second, complete independence of all counterterms on specific 
quark flavours implies that the 'equal flavour' and 'unequal flavour' ones are equal. E.g. 
[Xfo ir ]qgq'qi = [Xfo ir ] q gqq and so on. We have used this fact in writing the equation, hence, 
in (A.22), / and /' are not necessarily distinct flavours. Then using the flavour summation 
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rules, we can rewrite the summation over the unobserved indices k, t, i and r in the (m+2)- 
parton configurations into sums over indices kt and ir in the m-parton configuration. We 
have 

fc tj^k r^k,ti^k,t,r #(/kt)m #(firi i r kt) m 



(A.23) 



where the notation is the same as in eqs. (A. 9) and (A. 16). We assume that in general k, 
t, i and r are all distinguishable. Then we find 

#(9/)m+2 #(g;k ^ t) m+2 #(qf,r ^ k,t) m+2 #(g;i ^ k,t,r) m+2 _ . . 

(Trig + Z)\vrig + 1J , 







#(?/)m+2 #(£/; 7^ t) m +2 #(9; F 7^ fc) *)m+2 #(. 


g;i ^ M,r) m+2 _ 


#(?/)m #(5; ^ / fc*)m 




#0)m+2 #(5;^ / *)m+2 #(5;^ / k,t) m+2 #{g 


;i 7^ k,t,r) m+2 _ t 



(m g + 2)(m g + l), 
(m g + 2)(m g + 1) , 

■ (mp + + 1) , 

(mp + l)(mj, + 1) . 

#{g)m #{g; ir + kt) m 

(A.24) 

In case 3 we must remember that the counting is slightly different for the 'equal flavour' 
and 'unequal flavour' contributions even when the counterterms are the same. We have 

#(g/)m+2 #(?/; k / t) m+2 #{qf,r / k, t) m+2 #(<?/; i ^ k,t, r) m+2 _ 



#(g)m #(g; ir / kt) m 

#(q/)m+2 #(g;fc / ^)-m+2 #{qf,r / k,t) rn+2 #(q~f>;i ^ k,t,r) m+2 

#(?/)m #(5; / fc 0m. 
#(g)m+2 #(g; fe / *)m+2 #(gj*j r / M)m+2 / k,t,r) m+2 



#(g)m #(g; ir + kt). 

{m f + 2)(m/ + l)(m/ + 2)(m/ + 1) 



e.f. 



m 3 (m 9 - 1) 

#(qf)m+2 #(q~f,k / *)m+2 #(qf,r / k,t) m+2 #(qf,i / fc,t,r) m+2 



(A.25) 



#(s)m #(5; «> / fc*)r 

(m/ + l)(mj + l)(m// + l)(mj, + 1) 



u.f. 



mg(m g - 1) 

By exploiting the symmetry of this factor under permutations of indices, we trivially obtain 
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the rest of the cases as well. Finally, using eqs. (A.4)~(A.6), we have 



E ^EEE E 

{m+2} i m+2 i t k±t r^k,ti^k,t,r 



Egf-E E 



{m} 



kt ir kt 



ry(0) l(-) , fy(D] iW , ryW ll-J , fyW 1 
V^ktiriqgq'g " r l^fctirJoooo' ~r L^-fctirJooo'o ~r l^ktiri 



,(0) ,(...) 
*~ktir'qggq 



(0) ,(...) 
ktirigqq'g 



'(0) 

^ktir'gqgq' 



■ > fkt,q°fir,q 



+ {q^q) + (q' ^q) + (q^q,q ^ q) 



+ 



ry(°) i(-) 4- ry(°) l(-) 4_r> fry( 'iW , ry(°) i(-) 

l^ktirlqggg ^ktirlaqgg ' tf I i^ktirlqgq'q' V^k 



L ktir'qgq'q' 



+ 



fX -(0) iO..) , ry(0)i(...) \ 
y^ktirlgqq'q' l^ktirlgqq'q' ) 



+ 



ry(0) , r r (0) !(...) , fryfO) , ry(0)i(...) 

V^ktiriggqg i ktirigggq "T" "I ^KHtirVg'gg ^ i^ktiriq'q'qg 



, ry(0) , ry(0),(...) 

"T" l-^ktiriq'q'gq "I" l^ktiriq' q' gq 



S ht,9 S fir,q +{Q^ 



+ 



ryWlO-O 4-r< Ay(°)l(-) 4. ly(°) l(-) 4. fy(°) i(-) 

l^ktirigggg ~r" "f I ^ktiriggqq "i" t-^ktiriggqq ~r l^ktiriqqgg 



r l-^ktirlqqgg j ' l f I l^ktirl qqq'q' ktiAqqq 1 q> V^ktir\qqq< q' 



qqq'q' 



■ fy(0) \ 
I ktiri qqq'q' I 



J fkt,9°fir,g 



'qqq'q' 



(A.26) 



We remind the reader that eq. (A.26) was derived by using that the counterterms are 
independent of specific quark flavours (as the notation implies), and further that the 'equal 
flavour' and 'unequal flavour' subtraction terms are equal. Then the sums YJy and YJj/ in 
eq. (A. 22) may be performed, and we obtain the factors of rif and nj? as shown. 

Finally, we define the flavour summed counterterms as 



E 

{m+2} 



1 



'{m+2} 



EEE £ i^ 1 '--'- - - =v 



1 



V 



t k^tr^k,ti^k,t,r 



ktirifkftfifr Z-f Q r , Z.^ 



- |m| 

I" 1 } fct ir^fci 



E K 



'(0) \ 



(...) 



' Jktjir 



(A.27) 
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We obtain 



(0) 



> 



(...) 




ktir 



i ktiriqgq' g 



( X ktir) 





(A.28) 




(...) 



_ry(°)l(") 4_„ /ry(ll) l(-) i ryC) l(-) , r l(-) , fy(0) i(-) A 
~~ l^ktirlgggg t '*f ^rHtirJffSgg ~ r Y^ktiriggqq i^ktirlqqgg i^ktirlqqgg J 

, ^AytOjlW , r y (0) ,(...) r y (0) nC-) , f y(0) if...) \ 

"f ^L^fcjjrJgqg'g' ~r ktir 1 qqq> q> t l^ktiriqqq'q' l^ktirlqqq'q' J 



A. 2 Computing the flavour summed integrated counterterms 

Using eqs. (A.12), (A.13), (A.19), (A.20), (A.27) and (A.28), it is straightforward to com- 
pute all flavour summed integrated counterterms as presented in section 3.2, after taking 
account of the following points. 

• In eq. (3.6), most counterterms appear with some explicit overall factor, which must 
be included in the final result. E.g. for the collinear-triple collinear counterterm, this 
factor is 1/2. 

• In certain cases, the ordering of some flavour indices may be meaningless, due to 
a symmetry of the integrated counterterms. E.g. in the collinear-triple collinear 
case, the integrated counterterm is symmetric in the first two indices, [C kt C k °^ r ] qgq = 
[^kt^htr\ggg anc ^ so on - Hence, some terms that appear separately on the right hand 
sides of eqs. (A.20) and (A.28) may be equal. See also the second line of eq. (A.13), 
where the appropriate symmetry is already taken into account. 

• In particular cases, some terms that appear on the right hand sides of eqs. (A.13), 
(A.20) and (A.28) may be zero for certain flavour assignments. E.g. in the collinear- 
triple collinear case, the first two indices must correspond to a true singly-collinear 
limit, hence [C kt C$ r ] qqq - = and so on. 

• The 'equal flavour' and 'unequal flavour' counterterms in eq. (A.20) may actually be 
equal as e.g. in the collinear-triple collinear case, where [C kt C k °^ r ] qqq = [C kt C k °^ r ] qq i q i . 
Recall that this is already taken into account in eq. (A.28). 
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With these in mind, we easily find e.g. 

9 1 o o (A ' 29) 

= 2 ^kt^itrlggg + ™f [CfctCjjJ.]^ + 2nf [Ckt^itrlgqq > 

for the collinear-triple collinear flavour summed counterterms. The rest of the results in 
section 3.2 are obtained similarly. 

B Modified doubly-real subtraction terms 

We outline a simple modification to the NNLO subtraction scheme presented in refs. [1, 
2]. Parts of these modifications were presented previously: those relevant to the singly- 
unresolved approximate cross section dcr^'^ 1 appearing in eq. (1.3), and to the approxi- 
mate cross sections in eq. (1.4), were presented in ref. [35]. In this appendix we describe the 
modification of the iterated singly-unresolved approximate cross section do"^h^ 12 , which 
appears in eq. (1.3). 

Recall that the iterated singly-unresolved approximate cross section can be written 
symbolically as 

d<w' 2 Al2 = dcf> m [d P2 }A 12 \M { Z 2 | 2 , (B.l) 

where the iterated singly-unresolved approximation Ai2|Ai^j_ 2 | 2 is a sum of a number of 
different collinear-, soft-, and soft-collinear-type terms (see eqs. (3.1)-(3.4)). The precise 
definition of these terms involves the introduction of two momentum mappings 

Mn+i ^ {P}t ) , and {p} n+1 {p}W , (B.2) 

which are iterated in various combinations to produce appropriate mappings of m + 2 — >• m 
momenta. As discussed in section 2.3, all such mappings lead to an exact factorisation of 
the m + 2 particle phase space, symbolically written as 

#m+2({p};<3) = d(p m {{p}r^Q)[ d Pl,m][dpi,m+l] ■ (B.3) 

The exact form of the factorized phase spaces [dpi n ] (re = m, m + 1) is given in eqs. (2.16) 
and (2.17), but their only feature which is relevant presently is that they carry a dependence 
on the number of partons, n, of the form 

[dp£J] oc (1 - a vr f^~^ , (B.4) 

[dp&] * (1 - VrQ)^- 1 ^- 1 ■ (B.5) 

The subtraction terms, as originally defined in ref. [1] do not depend on the number of 
hard partons, thus the m-dependence of the factorized phase space measures is carried 
over to the integrated counterterms, where furthermore this dependence enters in a rather 
cumbersome way (see e.g. eqs. (A. 9) and (A. 10) of ref. [3]). 
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Thus, as in ref. [35], we reshuffle the m-dependence of the integrated counterterms into 
the subtraction terms themselves, where it appears in a very straightforward and harmless 
way, through factors of (1 — a) and/or (1 — y) raised to m-dependent powers. For easier 
reference, we gather the results in table 5, where together with the subtraction terms, we 
give the momentum mappings used to define the term(s) and the function which multiplies 
the original counterterm to produce the modified one. The / functions appearing in table 5 
are defined as 

f(zQ,z,p) =Q(z q -z)(1-z)-p. (B.6) 

The pattern of modifications is hopefully clear: if the factorized phase space ap- 
propriate to a given subtraction term carries m-dependence through factors of (1 — a) 
and/or (1 — y) respectively, it is multiplied by a factor /factors of f(ao, a, d(m, e)) and/or 
f(yo,y,d'(m,e)). We emphasise that the form of the exponents d(m,e) and d'(m,e), is 
actually fixed by the prescription in ref. [35] (see eqs. (3.2), (3.12) and (3.13) in particular) 
and the requirement that the modified subtraction terms should still correctly regularise 
all kinematic singularities. In fact, we must have 

d(m, e) = 2m(l — e) — 2do , and d'(m, e) = m(l — e) — d' , (B.7) 

where do and d' are the same constants which appear in eqs. (3.2), (3.12) and (3.13) of 
ref. [35], i.e. 

do = Dq + d±e , and d' = D' + d^e , (B.8) 

where Dq,D'q > 2 are integers, while d\,d\ are real. Also, the parameters Qo and yo must 
have the same values for all subtraction terms, including the singly-unresolved ones of 
ref. [35]. 

Finally, we note that the modifications introduced above do not spoil any of the can- 
cellations which take place among the original subtraction terms, hence the modified coun- 
terterms are still a correct regulator of all kinematic singularities. This is not particularly 
hard to check explicitly, and is actually a manifestation of the fact that the various mo- 
mentum mappings obey several conditions in soft and/or collinear limits. As these were 
discussed in ref. [1], we do not go into further details here. 

C Basic collinear, soft and soft-collinear functions 

Certain basic functions appear repeatedly during computations in this paper. They all arise 
as integrals of various simple factors over factorized collinear or soft phase space measures. 
Below we define and give explicit integral representations of these functions. Some, notably 
I, J and /C, have been considered previously in ref. [35], albeit in somewhat more general 
forms. For completeness, we present these here as well, although only in the special cases 
used in this paper. 
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Iterated collinear counterterms 


Subtraction term 


Momentum mapping 


Function 


< -fct L fctr 


M^{p} (h)C£ ^{p} ( " ? ' w) 


/(ao,a fct ,d(m,e)) 
x/(ao,afei 7: ,rf(m,e)) 




M^{p} (fci)Ci ?{p} (i? ' fct) 


/ (a , a& t , d(m, e)) 
x/(a ,a~_,d(m, e)) 


p pe(0,0) p p pq(0,0) 
L -fct U - , fet;r ' U kt^ir;kt^kt;r > 


M^{p} (fct)iz >{p} (? ' fet) 


/(ao,a fct ,d(m,e)) 
xf(yo,VfQ,d'(m,e)) 


n o(0,Q) r r o(0,0) 
u fci°fci I ^kt^rkt^kt 


W ^ {p} (fct) ^ {p} (S ' fct) 


f(a ,a kt ,d(m,e)) 
x/(2/o,2/fctQ> d '( m > e )) 


Iterated soft counterterms 


Subtraction term 


Momentum mapping 


Function 


c ^(0,0) <j ™j(0,0) 
°t L -irt > "tU>ir;t ! 

o /7 ~o(0,0) 
°i ( -i r tCO ir . t 




f(yo,ytQ,d'(m,e)) 
xf(a ,a~~,d(m,e)) 


c p o(0,0) c r9 ,,(0,0) 
«->tWrt rt 1 °t L0 ir;t rt ) 

c ^ r o c(0,0) C c(0,0) 


{p}^{p} (t) ^{pY ? ^ 


f(yo,ytQ,d'(m,e)) 
xf(yo,yrQ,d'(m,e)) 


Iterated soft-collinear counterterms 


Subtraction term 


Momentum mapping 


Function 


p /»(0,0) p q pq(0,0) 

f> ^o(0,0) 
i ( -irt L0 ir;i 


{ P } ^ {?}<*> C -14 {£}(->*) 


f(yo,ykQ,d'{m,e)) 
xf{a ,a~~,d(m,e)) 


r O o(0,0) ^ q p o(0,0) 
^rf->t c 'rt ! ^kt V-krt°rt > 

C rltC s (0i0) r <> <> (0 ' 0) 

C rt COlCfc r (O rt , Ufc t O t O rt 

C /ic c(0,0) 
CfctO t LO ir . t O rt 


{p> Ji>{p}Wi^ { pp.o 


f(yo,ytQ,d'(m,e)) 
xf(yo,y?Q,d'(m,e)) 



Table 5. The modified iterated singly- unresolved subtraction terms are obtained from the original 
counterterms (first column) by multiplication with an appropriate function (last column). Also 
shown are the momentum mappings used to define the subtraction terms (middle column). The 
/(zq, z,p) function is defined in eq. (B.6) while d(m, e) and d'(m, e) are defined in eq. (B.7). 



C.l Collinear functions 

When computing the integral of the azimuthally averaged Altarelli-Parisi splitting func- 
tions over the factorized collinear phase space, the following function arises: 

X i(^>Q' e ' a O' d o,fc) = l(yf rQ ,e,a ,d ,0,k,0, 1) = 

16-7T 2 f ( ir ) 1 ( C " 1 '* 
= —^-Q 2e / [ d H.m+i]— z r ,if(a ,a ir ,d(m,e)) , 
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which, as indicated, is simply the X function of [35], for the special values of parameters 
k = 5 = and gi = I (for their meaning see ref. [35]). The integral is over [dp^ +1 ], i.e. 
the factorized measure obtained when going from m + 2tom+l partons via the collinear 
mapping of ref. [1] , which can explicitly be written as 

[dpfi+M, P<r5 Q)\ = ^2-& (1 - a ir ) 2m ^~ l y frQ da ir dv ir 

x ap (a ir + (1 - a ir )y trQ )-* vp (1 - v ir y* ( C - 2 ) 

x @(a ir ) 0(1 - a ir ) Q(v ir ) 9(1 - v ir ) . 

The general collinear function was first computed analytically in ref. [36]. For convenience, 
we present the integral representation for the specific case, used in this paper: 



li(x;e,a ,d ,k) = x / da (1 - a) 2d °~ l a~ l ~ e [a + (1 
Jo 

x / dvv~ e (l — v)~ 



a + (1 — a)xv x k 
2a + (1 — a)x 



(C.3) 



where the factors 



i , /-i \ \ , a + (1 - a)xv 

a (a + (1 - a)x) and — — (C.4) 

2a + (1 — a)x 



correspond to the collinear pole yi r and momentum fraction z r .i, respectively (with x = x^_, 
a = ai r an d v = 1 — Vi r ). 

Among the iterated subtraction terms considered in this paper, we find two other basic 
integrals over the collinear phase space measure (C.2). One of these is the Lorentz tensor 

^(Pir,Q; e ,^0,d ) = ^—Q 2e I — z i,rZr,i "V/ ^ /(«0, OL ir , d(m, e)) 

(C.5) 

with kinematic dependence only on p^ r and Q^. The transverse momentum is defined to 
be orthogonal to both of these, Pi r -k±^ jr = Q-fe±,j,r = and contraction with g^ v replaces 
the fraction in the last factor with four. The most general Lorentz structure that obeys 
these conditions is 



1 - e 



nv _ PirQ" + Q^Pir , Q 2 ^fi 

9 ~ Pir-Q (p lr -Q) 2PlrPlr . 



X ^~Q 2e I [&P\%+\(PU Pir,Q)]— z r,iZi,rf(oiO, "in d(m,e)) , 



(C.6) 



where the integral in the second line is clearly just [Ii(xf r , e, ao, do, 1) — Ii(xf r , e, ao, do, 2)]. 
The other one is 

T 2 ix tr \ e, a , d ) = ^-Q 2e f [dp^ +1 ] — ^-f(a , a ir , d(m, e)) . (C.7) 
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The definition of Zi, r [1] implies yiQ = Zi <r (yiQ + y r g), and from the definition of p ir we 
have yiQ + y r Q = 2cti r + (1 — ot.i r )x<~. The momentum fraction can be expressed with the 
integration variables as in eq. (C.4), therefore, 

ViQ = a ir + (1 - a ir )xf r v ir , (C.8) 

hence 



rao 

l2(x;e,ao,do) = x / da (1 — a) 2d °~ 1 a~ e [a + (1 — a)x\ 
Jo 

x f di)u" e (l-i;)" £ [a + (l-QH" 2 . 



(C.9) 



C.2 Soft functions 

When integrating the eikonal factor Sik/(si r Sk r ) over the factorized soft phase space, we 
encounter the following function: 

16vr 2 f ,s s- h (c - 10) 

which, as shown, is simply the J function of ref. [35], for k = 0. The integral is over 
+1 ] , i.e. the factorized measure obtained when going from m + 2 to m + 1 partons via 
the soft mapping of ref. [1]. This measure can conveniently be written using energy and 
angle variables in the centre of mass frame, where 

(C.ll) 

p^ = E r (l, 'angles', sin sin sin 77, sin i? sin tp cos rj, sin $ cost/?, cos 1?) , 

where (here and below) the dots stand for vanishing components, while the notation 'angles' 
indicates the dependence of pr on the d— 4 angular variables that can be trivially integrated 
in all relevant cases. In terms of the scaled energy-like variable y r Q and the angular variables 
i9, f and rj, the two-particle phase space dfoipr, K; Q) is 



(C.12) 



dMPr,K; Q) = ^Ll_S e (- e ) 4 e dy y l ^5{y - y ?Q ) 

x d(cos ■&) d(cos <p) d(cos ry)(sin $)~ 2e (sin y?) _1_2e (sin rj)~ 2 ~ 2e . 
Often the integrand does not depend on all angles and we can integrate out 77, 

^ d(cos 77 )(sinr ? )- 2 - 2e = - — | _ ^ , (C.13) 

and 

ff r(ir 2 1) d(cos ^ )(sin v) ~ 1_2£ = 2 * • (ai4) 
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The J\ function of eq. (C.10) was first computed analytically in ref. [36]. We recall 
that it is conveniently evaluated in the frame eq. (C.ll), with the orientation fixed by 

p? = E 7 (l,...,l), p{J = ^(l J ...,sinx,cosx). (C15) 
The precise definitions of pf and p^ via the soft mapping [1] imply 

Sik = (1 - yrQ)S~l , S ir = S^ r , Sfcr = Sfc r , (C.16) 

and expressing all two-particle invariants with integration variables, we find [36] 

Ji(Y; e, yo, d' ) = -*Y ^ _ > n n (coa X (X), 1, 1) J dyy' 1 ' 2 ^! - y) d ° . (C.17) 



Above Y = y^-j, q, while 

cosx(^) = 1-2Y, (C.18) 
and the function £lji denotes the angular integral 

-l 



3 

' -1 



2e/„- ,_\-l-2e 



f2j7(cosx, Pi, fa) = / d(cost?) d(cos</?) (sini?) e (sin(p) 



(C.19) 



X (1 — Pi cosi?) ^ [1 — ^(sin x sin t? cos 93 — cos x cos 7?)] 



Presently we need the special case Pi = P2 = 1, which we call the 'massless' angular 
integral. The result of this angular integration is well-known [40] and is proportional 
to a hypergeometric function. Finally, using some hypergeometric identities and a one- 
dimensional integral representation of the hypergeometric function, we derive the following 
integral representation for Ji, to be used in this paper: 



Ji(Y;e,y ,d' ) = -Y-* v{ l_ 2e > J dy y- l ~ 2 %l - y) d ° J dt [l - (1 - Y)t 

(C.20) 



In some cases the eikonal factor involves three momenta as in 

is (ifc) ,(r) = Sil + SM . (C.21) 

+ SkrJSlr 

Then the soft integral 

J (lm HY (imQ ,P G% ye,y X) = 

167I " 2 ^2e fu (r) , Hi + Ski t < J I W ( C ' 22 ) 

= / [dPi,L+i(Pr,Q)\j———^—f(yo,yrQ,d{rn,e)) 

J e J \Si r -f- Sfc r JS/ r 

cannot be expressed with the soft function J any longer. In eq. (C.22) Pn^\ is the velocity 
of the momentum pf + p^ in the centre of mass frame. We evaluate J7^ lm ^ in the frame 
(C.ll), with orientation specified by 

pt + p^ = %£)(!,■■■,%£)), pf = E 7 (1,..., sin X , cos X ) ■ (C.23) 
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Using the precise definition of (j = i, k, I), we find 

su + s k i = (1 - VrQ){s^i + S£7) , Si r + s kr = s^ r + s^ r , s h . = s^ r . (C.24) 

Then we express all relevant two-particle invariants with integration variables and obtain 

. , r2(-\ _ f ) rvo 

J^(Y,P;e,yoA) = ~ W 2 J {1 _ ^ li(cos X (y, £), /M) J q dyy-^l-yf*, 

(C.25) 

where we used the soft phase space eq. (C.12) with rj integrated out (see eq. (C.13)). In 
eq. (C.25) Y = Y {U)%Q , p = p^, while 

cosx (Y,P) = ^^- (C.26) 

and the function Qji is defined in eq. (C.19) above. Now we need the special case Pi = (3 < 1 
and fa = 1) which we call the 'one-mass' angular integral. The evaluation of this integral 
will be discussed elsewhere, and here we simply indicate that the result is proportional to an 
Appell function of the first kind. Finally, using a one-dimensional integral representation 
of the Appell function, we obtain an integral representation for ^J^ m \ similar to eq. (C.17): 

2fi _ A rvo 



J^(Y, ft e, y , d> ) = -(2y)- 2 - 2 r r (1 ( ^ 2e £ ) ) jT dy y~^(l - y) d '» 



f dtr 1 - 
Jo 



l + P-(l + P-2Y)t l-p-(l-P-2Y)t . 

(C.27) 

Setting p = 1, we see that J'( lm \ Y , 1; e, yo, d' ) = Ji{Y; e, yo,d' Q ), as expected. 
C.3 Soft-collinear functions 

The following function arises when integrating the collinear limit of the eikonal factor, 
2zi tr /(si r z rt i), over the factorized soft phase space (recall that Zi ir /z r> i = SiQ/s r Q): 

Kt(e,y o ,d' o ) = K:(e,y o ,d' o ,0) = ^Q 2e [ [d^ +1 ]— ^/(y , VtQ, <t(m, e)) . (C.28) 

&e Jl Si r %r,i 

As indicated, this is just the K. function that was defined and computed in ref. [35], for 
k = 0. The integral is over [dp( m+ i\, i.e. the factorized measure obtained when going from 
m + 2 to m + 1 partons via the soft mapping. An integral representation for K, is easily 
derived in the frame of eqs. (C.ll) and (C.15) using 

Sir = s lr , s iQ = (1 - y r Q)Si Q + s~ ir , (C.29) 

which follow from the precise definition of p^ via the soft mapping. We find 

2(1 -y) 



IC 1 (e,y ,d' ) = 2^ dj/y- 2e (l-y) d »- 1 / d(costf) (sintf)^ 



o J-i 



1 + 



y(l — cos i?) 



(C.30) 
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which can also be written as 



rvo /■! 
K 1 (e,y ,d' ) = 2 dy y' 1 ' 2 ' {I - y) d ^ 1 dz z^(l - z)^(l - y + yz) , (C.31) 
io Jo 

where we made the substitution cos?? — > 1 — 2z. Notice that the soft-collinear function JC 
is independent of the kinematics. 

In some cases, the collinear limit of the eikonal factor involves three momenta as in 

2 1 - zt,ir 2 y iQ + y rQ ^ 



S(ir)t z t,ir S it + S rt ytQ 

Then the soft-collinear integral 

/C( lm )(% f) ;e,y ,do) = ^Q 2t [ [dp^ipf, Q)] — ±—^f{w, VtQ, d'(m, e)) 

(C.33) 

cannot be expressed with the soft-collinear function 1C any longer. In eq. (C.33), ^1?) * s 
the velocity of the momentum pf + ptf in the centre of mass frame. Using the definition 
of the mapped momenta [1], we have 

Sit + S rt = S~ H + Sr t , S iQ + S r Q = (1 - y t Q,) (s^ Q + S ? q) + S^ t + Sf t . (C.34) 

Then we evaluate /C^ lm ^ in the frame given in eqs. (C.ll) and (C.23) (with the trivial 
replacement k — > r). Expressing all two-particle invariants with integration variables, we 
find the following integral representation for the 'one-mass' soft-collinear integral: 



JC^(fJ;e,y ,d' ) = 2 2 * dy y-^{\ - y)^ 1 d(cos 0) (sintf)- 2e 



2(1 -y) 

1 + 



o J-i y(l -/3cosi?) 

(C.35) 

where f3 = Prq^y We make the substitution cos?? — > 1 — 2z to obtain the final form of the 
integral representation used in this paper: 

/C( lm )(/3; e ,,oX) = 2 rW 1 - 2 ^-^" 1 fd^l-^r 2 -;^ 

JO JO i — p + zpz 

(C.36) 

For /? = 1, we recover the soft-collinear integral, K.( lm \l; e, yo, d' ) = ICi(e,yo,d' ). 



D Integrating the collinear-type counterterms 

In this appendix, we discuss the integration of the collinear-type counterterms of section 4.1. 

(i) 

In particular, we define and give an explicit integral representation of all functions 
(t = l,...,9). 
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D.l Treatment of azimuthal correlations 

Collinear subtraction terms contain azimuthal correlations if the factorisation formula in 
the corresponding collinear limit involves gluon splitting kernels. For the two-parton split- 
ting ij — > i+j, these azimuthal correlations involve a transverse momentum k ± i ■ that we 
always chose to be orthogonal to the parent momentum p^. This condition is sufficient to 
prove that the integral of the spin-dependent and spin-averaged splitting kernels over the 
factorized phase space of the unresolved parton are the same [38]. Therefore, one can always 
substitute the spin-dependent splitting kernels P^J^Zj^kx^j'^) with their spin-averaged 

counterparts, {Pf°l.(zjj; k± t i t j; e)) = P^°J_(zj^;e), as done in section 4. 

In the case of strongly-ordered three-parton splittings, one cannot directly use the same 
argument. The splitting kernels in the integral of the collinear-triple collinear subtraction, 
C kt C^ t f^ depend on the transverse momentum kj_ kt in two ways. One is when the Lorentz 
index of the transverse momentum coincides with that of the parent gluon as in the explicit 
k± k t^± kt/h±kt terms in the gluon splitting kernels: 



9rqkqt V~«,Tj T,K> "jfc^f ) "f ,kt' -i-M,?' 1 

rr uu( Z r,H , Z kt,r , S rk ± , kit \ Z kt,r ^±,k,t^±,k,t 

2C A T R - g>* — 1 \-Zk,tZ t ,kT2 +4z k ,tZt,k ~ — p 

\ Z kt,r Z r,kt k ±,k,t S ktr / Z r,kt fc ±,fc,t 

4C A (1 - e)z 9mt ~ P£Uz k ,t,z t , k , fcx, M ;e) ^M^r M (m) 

r ,kt 



and 



^\Pg° g S}{z k ,t, zt,k,k± At , z £t? ,z ? r t , k £ ? ; e)\u) = 4C| 



Z r,kt , Z kt,r 



Z kt, r Z r .kt , 



1/ 



( Zk,t . Z tk \ U u 1_e S rk ±<kit Z kt,r k 1,k,t k ±,k,t 

X \ lTu TrA 9 z k,tz t ,k^T- 2 — - 2(1 - e)z kjt zt jk 

\Z t)k Z k .t J * K ±,k,t S ktr Z r,kt K ±,k,t 

- 4c A (i - ^ fii? ^(^v.^;0 ±,r / i ; ,tt • (D.2) 

±,r,fef 

This transverse momentum is not orthogonal to p ktr , as defined originally in ref. [1]. Never- 
theless, when integrating these subtraction terms, we can still substitute the spin-dependent 
splitting kernels with the spin-averaged ones, as we now show. 

Recall that the strongly-ordered three-parton splitting kernel appears in the collinear- 
triple collinear subtraction term in the form 

and the bra-ket expression above has the following precise meaning: 

(D.4) 
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where d^i (p ktr ,ni) and d vv < ( p ktr , n 2 ) are gluon polarisation tensors with (time- like) gauge 

fkftfr 



vectors, n\ and n 2 - Hence, the object to integrate is not simply (fi\P s f ° f y \u), but rather 



the contraction 

d^'{p ktr^ n l){AP}'^!;fM)dv'v{P ktr^2) , (D.5) 

that is clearly orthogonal to p ktr (because of the presence of the polarisation tensors). 
Thus, by the usual arguments, the azimuthal correlations present in C kt C^^ vanish after 
integration over the phase space of the unresolved parton. However, we must still be careful 
to compute the average over the polarisations correctly. When k±- p 7^ 0, we have 

a (~ X^'w (~ \ 1 (v k± -P \fi> k± -P \j (r>r\ 
d^'{p,ni) ^ 2 d v , v (p,n 2 ) = -rj- I k ±ill - — I kx )V - — n 2v ]+•••, (D.6) 

K± Kj_ \ p -Tlx J \ p •Tl2 J 

where the dots stand for terms proportional to p^ or p v which vanish after contraction 
with the matrix element, by gauge invariance. Thus we find (n is a further time-like gauge 
vector) 

/, ,~ ^± k ± , ,~ A 1 , / ,d tllxl (p,n 1 )k^kfd v > u {p,n 2 ) 

( V(P- n i) ,2 d v/v {p,n 2 ) ) = Try- r<v(p, n)- 



1 



2(1 - e) k\ 



k±-p fci-p, (k±-p) 2 

— fcj_-ni - — k±-n 2 + rTTc rni-ri2 

p-ni p-n 2 {p-n\){p-n 2 ) 



(D.7) 



Eq. (D.7) shows that the advantage of having k±-p = actually lies in the trivial azimuthal 
averaging. However, this can also be arranged if k±-p ^ 0. For example, choosing gauge 
vectors and such that + oc k^_ (and of course n\ = n 2 = 0), we find that the 
above average just reduces to — 1/[2(1 — e)], which is the usual result. Since k\ < 0, such 
nonzero and always exist. Indeed, in any Lorentz frame we have the parametrisation 
k^_ = K(l,j3v), where v 2 = 1 and f3 2 / 1. Then setting e.g. 

nZ = ±K(l + P)(l,H), and n£ = \k{1 - /3)(1, -if) , (D.8) 

we clearly have and nonzero, n\ = n 2 = 0, and + = k'f , as required. We make 
essentially such a choice in our numerical code. 

A shorter, though less transparent proof is to observe that if we change this single 
troublesome term to 



where 



then obviously kj_,k,fPktr = 0> an d an the usual arguments apply. What is not immediately 
obvious, is that this modification does not ruin any of the delicate cancellations between 
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the various counterterms in any IR limit, and hence is allowed. This can be established as 
follows: in the centre of mass frame of Q, the replacement of eq. (D.9) simply amounts to 

kj_ k ,t = -K(l,(3v) -> fcj_ iM = ~K(l + A, (3v) , (D.ll) 
with some A. Then, continuing to choose the gauge vectors as in eq. (D.8), we have 



< = \k(\ + P)(l,v) -+ nl = X -K(\ + A + f3)(l,v) , 
n% = Ik(1 - /?)(!, -v) -> n» = \k{1 + A - /3)(1, v) , 



(D.12) 



i.e. only the normalisation of the gauge vectors is changed by the replacement in eq. (D.9). 
However, this implies that eq. (D.6) is actually unchanged. Indeed, recalling that = 

. i* - .K' + ^')K' + <) 

(ni + n 2 ) 

(D.13) 



[p,n 2 ) = a^'{P,ni) — ± — , „ N2 

(p,n 2 ) 



it' 

2(ni-n 2 ) 

and the last expression is clearly seen to be invariant under (independent) rescalings of n\ 
and n 2 . Hence, the replacement in eq. (D.9) is completely harmless. (However, notice that 
the proof requires the specific choice of gauge vectors as in eq. (D.8).) 

The other occurrence of k± j k,t m the strongly-ordered kernels is in the ratio — -s — — , 

also present in the quark splitting kernels. Examining the explicit forms of the strongly- 
ordered splitting kernels, we find that this ratio always appears in the form 

PsKl^(p kt ,Q), (D.14) 

ktr 

where the integral is defined in eq. (C.5) and computed in eq. (C.6). Contracting the 
latter with PrPr/ s £tr J we obtain 

2 f Q 2s £tr s tq\ 16vr 2 2e f (fct ) , 1 



\ ktQ kt Q / 

Observing that 



Q 2€ f[dp^ +1 ] — zt,kZk,tf(ao,(Xir,d(m,e)). (D.15) 

J ' Skt 



^ktQ kt,r 

we find that when integrating the strongly-ordered splitting kernels over the factorized 
phase space, the integrals of 

1 2 6 k? ?k± ' k ~ and (D.17) 

-L,k,t^ktf ktQ kt,r 

are equal, so we can substitute the former with the latter, which we implement in the next 
subsection, where we give the spin-averaged splitting kernels explicitly. 
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D.2 Explicit forms of the spin-averaged splitting kernels 

For the sake of completeness we present the explicit expressions for the spin-averaged 
splitting kernels used in the integrals of the subtraction terms. 

The azimuthally averaged Altarelli-Parisi splitting kernels are well known: 



2C A 



1 1 

_+ i 

z 1 — z 



2 + z 



p£\z;e) = T R 



P$(z;e) = C F 



--2 + (l-e)z 



(D.18) 
(D.19) 
(D.20) 



In our convention the ordering of the labels on the splitting- kernels is usually meaningless, 
but in eq. (D.20) z refers to the momentum fraction of the second label. In other words 
Pgq\z; e) = Pqg\l — z; e). The other two cases are symmetric with respect to z ■<->• 1 — z. 

In the strongly-ordered kernels Pf k °f}f^ the ordering matters, too. As a result, the same 
triple-parton splitting function may have different strongly-ordered limits, which can be 
distinguished by the momentum labels in the kernel, once the ordering of the limits is fixed 
by the momentum mapping, ktr — > kt + r — > (k + t) + r in our convention. We always 



choose z = ztk and 3" 



■M 



as independent variables. For quark splitting we have 



Pqgg ^ ^ ( z i z ) e ) ~~ ^qg ( Z ) e )-^gg ( z i e ) ; 

P^-(°)(z, z,R;e) = P$(z)P$(l - *;e) - C A C F z(l ~ R, 



P s ;J (0) (z, z,R-e) = P$(z;e)P$(l -z;e)- C A C F z(l - z)b^R 



q'q'q v~: ~ ) — ) -/ - qq 

while for gluon splitting we find 

P^(z,z,R) = P$(z)P$(z) - C\z(l - z)(b g ffR, 

P s q£ i0 \z, z,R;e) = P^\z;e)P^(z) - C\z(l - z^b^R, 
P s gq % (0) (*, z ■ e) = P$ (1 - z; e)P$\z; e) , 

where the constants b*\ are given in eq. (4.3). Eqs. (D.21)-(D.26) can also be written in 
a unified form, 



(D.21) 
(D.22) 
(D.23) 

(D.24) 
(D.25) 
(D.26) 



,(0) 



,(0) 



fkfM'^ P fktfr 



We see that the second term is present only if the three-parton splitting involves a two- 
parton sub-splitting with parent gluon. 

D.3 Collinear-triple collinear counterterm 

The collinear-triple collinear counterterm involves two successive collinear mappings, which 
leads to exact phase space factorisation in the iterated form 

d<p m+2 ({p};Q) = dcj> m ({p}^ t? ^;Q)[dp^\p r , p% ; Q)][dp^ +1 ( Pk , Pkt ;Q)] . (D.28) 
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The one-parton factorized phase spaces [dp^!^ l+1 (pf.,p kt ;Q)] and [dp^^\p r , Pg t ^; Q)] 
are given explicitly by eq. (C.2) after appropriate changes in labelling, including the re- 
placement m — > m — 1 in the second case. The Altarelli-Parisi splitting functions and the 
factor z(l — z) in eq. (D.27) can be expressed as linear combinations of powers of momen- 
tum fractions. Consequently, the integral over the factorized phase space [dp^ +1 ] in the 
integrated collinear-triple collinear counterterm is written in terms of collinear functions 
li(xg t , e, ao, do; k) of eq. (C.l). In order to compute the subsequent integrals over [dpj*^], 



#\ t?\ zf } = J* / [dpg;']^ (i - Ct) >, j 



tr 



ktr \ y£tQ) (D.29) 

x f (&o, ag_ ? , d(m, e))li(xg_, e, a , d ; k) , 

the variable X£ t = y£ t q needs to be expressed in terms of p^ tr instead of p£ t (see eq. (C.8) 
with proper changes in labelling), 

X kt = VktQ = a ktr + (! - a ktr ) X Wr V kt r ■ ( D - 30 ) 

Then using the abbreviations a = ol^, v = vr.^ (the integration variable corresponding 
to 1 — z~ £ t ) and x = x-rr , the integral representations (C.3) and (C.9), we find that the 
integrated collinear-triple collinear counterterm can be expressed using the following three 
types of integrals: 

l^\x;e,a ,do;k,l) = x / da(l - a) M °~ 3+2e a _1_e (a + (1 - a)x)~ l - e 

Jo 

f 1 c , x c f ot + (1 — a)xv \ 1 ,„ 
x / dvv- e (l -v)~ e - v '—) (D.31) 



Jo \2a + (1 — a)x j 

x T\{a + (1 — a)xv, e, ao, do; k) , k, I = — 1, 0, 1, 2 , 



lP(x;e,a ,d ;k,l) =x / da(l - a) M °- 3+2e a _1_e (a + (1 - a)xy l ~ e 

Jo 



f rl x r ( a + (1 - a)x(l - v)\ l 

x / dvv~ e (l-v)~ e [ ' \ - D.32 

Jo V 2a + (l-a)x J V ' 

x X\{a + (1 — a)xv, e, «o ; do; k) , k, I = —1,0, 1, 2 , 



and 



xl 3) (x;e,ao,do;/c) = x I " da(l - a) M °~ 3+2£ a- £ (a + (1 - a)z) _ 



oil 





x / di)t)- £ (l-!))- £ (a + (l-a)ra)" 2 (D.33) 
Jo 

x Xi(q + (1 — a) xv, e, ao, do; k) , A; = 1,2. 
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We can use the relations 



and 



xj?\x;e,a ,d ;k,0) = X^\x;e,a ,d ;k,0) , (D.34) 
xf?\x; e, a ,d ; k, 1) = T$\x; e, a , d ; k, 0) - X^' \x; e, a , do; 1) , (D.35) 

xf\x; e, a ,d ; k, 2) = X^\x; e, a , do; ft, 0) - 2X^\x; e, a ,d ; k, 1) 

+ X^\x;e, a ,d ; k,2) , (D.36) 

to reduce the explicit computation of the X c integral to the case I = — 1 . 

In terms of the functions X^ (i = 1, 2 and 3) we find the result given in eq. (4.2). 

D.4 Collinear-double collinear counterterm 

The collinear-double collinear counterterm also involves two successive collinear mappings, 
which leads to exact phase space factorisation in an iterated form similar to that in 
eq. (D.28). The integral over the factorized phase space measure [dpf^+J leads to the 

same collinear integrals as in eq. (C.l). Then the necessary integrals over [dp^ 1 ^] are 



Xf = l ^~Q 2 ' [ [d^]-UU, f(a ,a 7? ,d(m,e)) 

€ J j r 

x Ti(xfo, e,a ,d ;k) . 



(D.37) 



Again, xc. needs to be expressed in terms of p^ t instead of p^ t , 



ifi = 2M = t|M = (1 _ aif)ls , (D .38, 



Setting a = ct^~, v = v^~, x = x-rr, y = xj^, and using the v 1 — v symmetry of the 
integration measure, we find that the integrated collinear-double collinear counterterm can 
be expressed as a linear combination of the integrals 



1-/11 f a ° 
X^>{x,y;e,a ,d ;k,l) =y / da(l - a) 2 *" 3 ^"^^ + (1 - a)?/]" 1 

Jo 



x / dvv~ e (l — v)~ 







a+(l-a)yvy (D.39) 
2a + (1 — a)y / 



x Xi((l — a)x, e, ao, do; ft) , ft, Z = — 1, 0, 1, 2 . 
In terms of the functions X^\x, y; e, a^, d^; ft, Z) we find the result given in eq. (4.7). 
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D.5 Collinear-soft-collinear-type counterterms 

The collinear-soft-collinear-type counterterms in eqns. (4.8), (4.12) and (4.10), involve a 
collinear mapping followed by a soft mapping of the phase space, which leads to an exact 
factorisation of the original m + 2-particle phase space in the form 

d0 m+2 (M) = d^ m ({p}( ? ' fc '))[dp^ (D.40) 

where the factorized phase space measures [dp^ +1 ] and [dp^] are given in eqs. (C.2) 
and (C.12) after appropriate changes in labelling, including the replacement m — > m — 1 
in the second case. The integral over [dp^2 +1 ] gives the same collinear function as in 
eq. (C.l). In order to compute the subsequent integrals over the measure [dp^], 



kt,r 



S 7r Z r,i S ktr Z r,kt J (D.41) 

x f(yo,yrQ,d'{m,e))li(x£ t ,e,a ,d ;k) , 

we have to express the invariants of the dependent momenta (with hat) with those of the 
independent ones (with tilde): 



s^r = C 1 -y?Q) s jT > s fc? = s fcf for k = i,j,kt,i, 

S kQ = ( l -yrQ) S ~ k Q + S ~kri for k=7,kt, 



(D.42) 



which also implies 
Furthermore, 



x 



kt 



0--V?Q) x kt+Vkt?- (°- 43 ) 



z i,r via ( l -yrQ)y~ i n + y~ i - 1 



(D.44) 



VrQ y?Q 

with a similar expression for zc. ~ / z~ r-. 

To write explicit integral representations of X^ 1 (i = 5, 6 and 7), we choose the specific 
Lorentz frame of eq. (C.ll) with a different orientation for each function. 

Integrated collinear-soft-collinear counterterm. Here and in the following, we will 
use the partial fraction identity below to disentangle the singularities associated with the 
factors of l/s^j~ and l/sy~, appearing in the eikonal factor (first term in the braces in 
eq. (D.41)): 

s 3l _i-y?q y~fi _ 1 -y?Q 4y __ y~jQ y ~lQ 



_i ( D -45) 

1 - vtq iY ( yjQ i y~iQ \ ( 2 yj? , 2 yiA 
Q 2 ll ' Q Vyi? 2 yir)\yiQ viq) ' 

This is useful when computing the integral via iterated sector decomposition, while the 
original form is better suited to derive the Mellin-Barnes representation. Since 'undoing' 
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the partial fractioning is trivial, below we will show the more elaborate form of the integrals, 
which are directly suited to treatment with sector decomposition. 

The convenient frame for integrating the first term in eq. (D.45) is 
P*i = E l0-, ■ ■ ■ > !) i Pi = Eji 1 ,--- ,sinx7,cosx7) , 

(D.46) 

P kt = E kt (!> • • • > sin ^ kt sin X kt , cos <^ sin x Tt , cos x S ) , 

while for the second term we choose a frame where j and / are interchanged as compared to 
eq. (D.46). In terms of the scaled energy-like variable y-fQ an d the angular variables #, <p 
and 77, the two-particle phase space d^2(Pr> ^! Q) * s give 11 by eq. (C.12). The two-particle 
invariants y~~, ^Xt f nave to be expressed in terms of the integration variables, i.e. 



-f-^=y fQ (l-costf), (D.47) 
y jQ 

— = 2/f q(1 — sinxy sin # cos </? — cosxy cos 1?) , (D.48) 



Via? = ^TtQVrQi 1 ~ sin <^ sin X^ sin sin <^ cos 77 

z (D.49) 

— cos </>^ sin x^; sin $ cos — cos x^: cos 1?) . 



Furthermore, writing out the definition, eq. (2.23), of all Yjj. q's in the specific Lorentz 
frames, we easily find that the fixed angles can be expressed with invariants as 



cos xy = cos xO^Jq) , cos XTt = cosxO^q) , (D.50) 

cos^ =co S( f>(Y llQ ,Y 1Tt>Q ,Y 1TtQ ), (D.51) 

with 



cosx(^) = l-2F, sm X (Y) = 2y/Y(l-Y), 

Y 1 +Y 2 -Y 3 - 2Y X Y 2 ( D - 52 ) 



COS 0(^2,^3 



2 y /Y 1 (l-Y 1 )Y 2 (l-Y 2 ) 



Using (C.12) and the expressions for the ratios of kinematic invariants in eqs. (D.47)-(D.51), 
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(D.53) 



(5) 

we find the following explicit expression for T£>\ 

, , /_2 2e e \ rvo , 

4 5, (x,yi,y 2 ,y 3 ;e, oo.do,^,^; A;) = - ( J 4y i J dyy-^l-y)^'^ 

x J d(cos i?) d(cos <p) d(cos 77) (sin $)~ 2e 
x (sin^)~ 1 ^ 2e (sinr/)^ 2 " 2£ (l -cosii)^ 1 
x ^2 — cos 1? — sinx(Yi) sin d cos (p — cos xO'i) cos "d\ 
x Jif[(l -y)x + y-^ t9 (y, a;, yi,y 2 ,y 3 ,^, 93)]; 6,00,^0;^ 
+ -y)x + y^ ?f (y,a;,yi,y3,y2,«?,^)];e,ao,do;A;^ , 

which we need for fc = —1, 0, 1, 2. Above, x = x-rr, Y\ = y~7™ Yo = Y~~ „, and 

kt ' 3 I ,Q j kt ,Q' 

y3 = yy q- The second term in the squared brackets in eq. (D.53) corresponds to the 
second term in the partial fractions of eq. (D.45). The integral representations of the two 

terms are formally identical, only the kinematic variables Y2 and Y3 are interchanged. 

(5) 

In terms of the functions X c we find the result given in eq. (4.9). In writing eq. (D.53), 
we have tacitly assumed that p 1 ^, pf and p^ t are all distinct (massless) momenta, whose 

kinematics is furthermore unconstrained. For processes involving two or three hard final 

(5) 

state partons, this is not the case, so the integral Iq with full kinematic dependence, as 
written above, first appears in computing NNLO corrections to processes with at least four 
hard final state partons. 

When there are only three hard partons in the final state, the kinematics of the event is 
constrained because momentum conservation forces the final state momenta to be coplanar. 
Thus, in eq. (D.46) we have sin^^r = 0, and the parametrisation of pj, p t and p kt 
simplifies accordingly: 



Pkt = E kt^ ■ ■ ■ >- sin XTt> cos Xkt) > 



(D.54) 



where we choose cos^^j = —1, so that we may assume sinxy and sinx^ to be non- 
negative. As a result of the constrained kinematics, we can first of all perform the cos 77 
integration in eq. (C.12) using eq. (C.13). Second, out of the four kinematic invariants (xj-r, 
Y-jj q, Y-.-j^ q, and Y-^ g) of the general case in eq. (D.53), only two are independent. 
(Momentum conservation implies three constraints among the five variables E~ , Ej , E-r: , 
cosxy and cosx^-) However, it is convenient to leave the formal dependence on all four 
variables, with the constraints 



Y l kt,Q Y j l,Q + Y j kt,Q 2Y j l,Q Y j kt,Q + 2 \/ Y j 1,q0- Yj l,Ql^j kt ,q0- Y j kt ,Q> 



(D.55) 
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and 



Y~j „(1 — Y-y 



kt 



3 kt 



' * 3 kt,Q\J Y j !,q( 1 Z,q) 



(D.56) 



The physical region for Y~j „ and y~n ^ is < Y~j „, F~-rr o < 1 and Y"~t „+Y~7t ^> > 1. 

J 4 ,Q 3 kt ,Q — 3 I ,Q J j kt ,Q — 3 I ,Q 3 kt ,Q — 

The two-particle invariant becomes independent of <^tt, 



% = 2^^ fQ(1 + Sin ^ Sm,5cOSV " COSXsC08,?) 



(D.57) 



therefore, the dependence on Yjy g, Y^YtQ an< ^ ^IktQ en ^ ers on ^y through the angles 

cosxy = cosx(yj T)Q ) , cosx^ = ^xO^q) (D.58) 

for the first term in the partial fraction and the angles 

cosxj = cos X (Y-j l Q ) , cosx^ = cosx(Y m Q ) (D.59) 

in the second one (with eq. (D.52) for cosx(y)). Integrating out cos 77 as in eq. (C.13), the 
integral (D.53) simplifies to 



(5) 

X <Ui( x > y i' y 2,^3;e,ao,do,yo,d ; fc ) = ~ 4y i 



r 2 (i- e ) 

2vrr(l - 2c) 



I/O 



dyy- l ~ 2e {l-y) d '^ l+e 



x / d(cos^)d(cos^)(sin??)- 2£ (sin^)- 1 - 2e (l-cos^)- 1 



x ^2 — cos i? — sinx(Yi) sin ■& cos (p — cos x(Xi) cos $ 
x XiM(l - y)x + y^ ? (y,x,Y 2 ,'&,cp)];e,ao,do;kJ 



(D.60) 



(/c = —1, 0, 1, 2) and it replaces the function 1^ in eq. (4.9). 

Further simplifications emerge if the three labels j, I and kt are not all distinct and 
e.g. I = kt, which is the only case relevant for processes with only two hard final state 
partons. (Note that j ^ I, so up to j I interchange, this is the only option.) Then 



Y — 

I kt( 



Yr-r 



I u 



x-n and Y~~r 



kt 



3 I 



eq. (D.46) simplifies to: 



and the integral (D.53) depends only on two kinematic variables, 
Q = fci Q ' P arame t r i sa ti° n °f the two hard momenta p j and p kt in 



p£ = E ? (l,...,l), 



Pkt = E kti l i ■ ■ ■ > sin X^,cosx 



■kt- 



SO 



'kt' 



2 y ktQyrQ( l ~ sin Xs sin ^ cos v? - cosx^ costf) 



(D.61) 



(D.62) 



-67- 



and the dependence on Yj^q enters through the angle 

cos XTt =cos X (Y irtQ ) (D.63) 
for the first term in the partial fraction and through 

cosxj = cos x( Y ]te )Q ) (D.64) 

in the second one (with eq. (D.52) for cosx(^O)- Integrating out cos 77 as in eq. (C.13), the 
integral (D.53) reduces to 

4%.(x,yi;6,ao,do,W),do;A:) = -4Yi ^ _ ^ J dy y"^ 2e (l - 

x J d(cos^)d(cos^)(sin??)- 2£ (sin^)- 1 - 2e (l-cos^)~ 1 

x ^2 — cos 1? — sinx(Yi) sin 1? cos cp — cosx(^i) cosfi^ 
x Xi ([(1 - y)x + y^ ? (y, x, Yi, 1?, e, a Q ,d ; k^j 

+ ([(! - y) x + Vktrivi x ' °' ^ ¥>)]; e, "0, do; kj , 

(D.65) 

(fc = —1, 0, 1, 2) and it replaces the function X c in eq. (4.9). 

Rest of the integrated collinear-soft-collinear-type terms. The convenient frame 
for integrating the second term in the braces in eq. (D.41) is 

Pi =E 1 (1,...,1), p^ = ^(l,...,smx^,cosx^). (D.66) 
The two-particle invariants are expressed with the integration variables as 

^ = yfQ (l-costf), (D.67) 

yktr = 2 y ktQ y?Q ^ - 5in Xg sintf cosv? - cosx^ COS 1?) , (D.68) 

where 

co SX5 =cosx(% Q ), (D.69) 

with eq. (D.52) for cosx(^)- Then, using eq. (D.44) the integral in eq. (4.12) can be ex- 
pressed clS £1 linear combination of the integrals 

ik\ r 2 i'i — a ry° 

Xf (x, Y;e, a , do, 2/0, d' ; k) = ^ _ ^ J dy y~ 2e (l - yfo'^ 

x / d(cos^)d(cosy) (sin^-^Csiny)- 1 - 26 2 I_^l±^) 
J -1 y(l-costf) 



Xli[[(l - y)x + y^ ? (y,x,Y,$,<p)];e,a ,d ;k) , k= -1,0, 1,2, 

(D.70) 
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as in eq. (4.13). 

Finally, the integral of the third term in the braces in eq. (D.41) is obtained by setting 
i = kt in the previous case, which implies Yr ^ _ — > Yr-r g = 0, so the integration over (p 
can be evaluated using eq. (C.14) and the integral in eq. (4.10) can be expressed as a linear 
combination of the integrals 

4 6 \x,0;e,a ,d ,yo,d' ;k) = 2 2 * / dy y~ 2 \l - y) d '»~ 2+e 

Jo 

x [ 1 d(cos tf) (sin *)-* 2 ~f l + C °f } (D.71) 
J-i y(l-costf) 

x XiT[(l - y)x + y^ ? (y,x,0,'&,ip)];e,a o ,do;k] , 

(k = -1, 0, 1, 2) as in eq. (4.11). 

D.6 Integrated collinear-double soft-type counterterms 

The collinear-double soft counterterm is defined by an iterated application of a collinear 
and a soft momentum mapping, which results in an exact factorisation of the original 
(m + 2)-particle phase space very similar to that in eq. (D.40). The difference is that in the 
present case the soft measure involves the momentum p^ t instead of p!^. The integrand of 
the collinear integral is the spin-dependent splitting kernel of gluon splitting, with Lorentz 
structure 

1 / i i m \ uiy (0) ( \ ,(0) 1-e ik ±,k,t k ±,k,t 

tHw p f k ft( z k,t, z t ,k,k±, k ,t; e)M = -r a M z t,k) - bV ft —r-z k ^ k — -2 , 

° A z K ±,k,t 

(D.72) 

where 



and in?/ t is defined in eq. (4.3), hence, the integral over the phase space measure [dp^ +1 ] 
involves both collinear functions I and X^ v '. Introducing the abbreviation 



X hh = ~^-Q 2e [ [ d PiS+i(P*> PktiQ)]— f(ao,akt,d(m,e)) 

be J Skt (D.74) 
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and using the result in eq. (C.6), we find that 

X hh = sT^f-CP [ [dp^ +1 ]^-f(a ,a kt ,d(m,e)) 

x { a hfM^) + b^ ft Zk,tz t ,k) 



■2 



, w i r lit -v ■ -* i- ret _ ^ -j/ 



(D.75) 



16tt' 



Q 2e / [ d Pi^2+i] — z k jZt, k f(a ,a kt ,d(m,e)) 
J s k t 



S, 

where in the first integral we recognise the spin-averaged splitting kernel, 

a hfi ^ + b hft Zt ^ 1 ~ z ^ = P hft ^ e ) • ( D - 76 ) 
Therefore, Zf k f t can be expressed as a linear combination of collinear functions I, with 
Lorentz structure exhibited in eq. (D.75). After contraction with S~^(kt) (see eq. (4.14)), 
we obtain 

l -S^(kt)X hft , y = - l -S v (kt) l ^Q^ |[d^ +1 ]i-P^ t (^; e )/(a ,a fet ,d(m, e )) 
. (0) / S 3Q s Tq 2Q 2 \ 



+ b 



fkft \ sr^sc^ sr,7Sc^ si 



ktj ktQ ktl ktQ ktQ , 



16vr 2 



Q 2e j [ d Pi^m+i] — ZtA 1 - z t,k)f(ao,u k t,d(m,e)) . 

J &kt 

(D.77) 



In the integral of the complete collinear-double soft subtraction, eq. (D.77) is multiplied 
with terms that are symmetric with respect to the interchange j I, and summed over 
both j and I (cf. eqs. (3.2) and (3.6)). Therefore, the integrals of the terms 

*l Q and SlQ (D.78) 

S ki3 S ktQ S ktl S ktQ 

in the parenthesis give identical contributions, and it is sufficient to evaluate three types 
of integrals: 

ler \,2> lu («>. 1 / "JT 2s Je 2l ? 2 l 



i«",4 8 ',i™} = ^</[d P g] 



S ktr [ S 3kt S lkt S ktJ S ktQ S ktQ ) (D.79) 

x f(yo,y ktQ ,d'(m,e))l(xg t ;e,a ,d ;k) . 

To compute the integrals in eq. (D.79), first we have to express the invariants of the 
dependent momenta with those of independent ones: 

sjj = (1 - VktQ^ s ]l ' s Jkt = s jkt> s lkt = s lkf 

(D.80) 



S jQ ~ ^ ^Mq) S jQ^ S jkf S IQ~^ y£tQ^ S lQ + S l 



kt 
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Then we can proceed as usual by writing the factorized phase space [dp^ +1 ] explicitly. 
We choose a Lorentz frame defined by 

^ = ^(1,...,!), pf = ^(l,...,sinx T ,cosx T ), (D.81) 
where the two-particle invariants are expressed in terms of the integration variables as 

1^ = y£tQ^ 1 ~ cos $) ' = yktQ^ 1 ~ sin xy sin ^ cos - cosxycosi?) , (D.82) 

y jQ y lQ 

with 

cos x~i = cos xQjT.q) ' (D.83) 

and eq. (D.52) for cosx(Y). 

For the eikonal factor (first term in the braces in eqs. (D.84)-(D.86)) we use the partial 
fraction identity (D.45) (with the substitution r — > kt), and the j I symmetry to 
integrate only one term. Thus we are left with the following three integrals: 

■p2 / 1 .A pyo 

lP(Y;e,ao,d ;k) = -4Y ^- J- j dyy" 1 " 2 ^ - yf^ 

x [ d(cos^)d(cos^)(sini9)- 2e (sin( / 9)^ 1 " 2£ (D.84) 

J_l 1 — COS 77 

x ^2 — cos i? — sin xOO sm $ cos <p — cos x(Y) cos i?) 
x Xi(y;e,Q ,d ; A;) , fe= -1,0,1,2 

(here Y corresponds to Y~. j g), 

4 8) (e,a ,d ;k) = 2 2 * f ° &yy~ l ^{l - y) d '^- 
Jo 

x f 1 d(cos^) (sinfl)- 2 * 2 "^ 1 *^ (D.85) 
7_i 1 — COS V 

x lx(y; e, a , d ; k) , k = —1,0, 1,2, 

and 

X ( l 9) (e,ao,rfo;fc) = -2 r( ^_~ 2e ; ) y o dy ^^(l - y) d °~ 2+e (D.86) 

xIi(j/;e,ao,d ^), k = 1,2. 

(The cases A: = —1,0 in eq. (D.85) are needed for eq. (4.16), see next paragraph.) The final 
result is presented in eq. (4.15). 



-l 
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For the integrated collinear-triple collinear-double soft subtraction in eq. (4.16) the 
phase space factorisation is the same as for the integrated collinear-double soft subtraction, 
and we need to evaluate the integrals 



Q 2e [« 



(fct)i j 2 z r,kt 2Q 



2 



' 2 ( 

\tr Z kt,r S ktQ J (D.87) 



x f(yo,y£ t Q,d'(?n,e))l 1 (x^ t ;e,ao,d ;k) . 

Recalling the definition of the momentum fractions, in the first term we recognise twice the 
second one in eq. (D.79) (after replacing j with r), while the second term is equal to the 
last one in eq. (D.79). Thus we do not have to compute any new integrals, we can express 
this integrated counterterm as a linear combination of the integrals defined in eqs. (D.85) 
and (D.86). The final result is given in eq. (4.17). 

E Integrating the soft-type terms 

In this appendix, we discuss the integration of the soft-type counterterms of section 4.2. 

(i) 

In particular, we define and give an explicit integral representation of all Tg functions 
(i = l,...,12). 

E.l Soft-triple collinear-type counterterms 

There are three integrated counterterms that involve a soft mapping of the momenta, 
followed by a collinear one, which leads to an exact factorisation of the original m + 2- 
particle phase space in the form 

d(f> m+2 ({p}) = d0 m ({p}(f ? '*))[dp£ ) (p r , PirlQWp^falQ)] • (E.l) 

The factorized phase space measures [dp^ m+l (pf, Q)] and [dpj 1 ^ (p r , p j r ; Q)] are given 
in eqs. (C.12) and (C.2), with appropriate changes in labelling and the replacement of 
m — > m — 1 in the second case. 



Integrated soft-triple collinear counterterm. We begin by noting that the soft func- 
tions pj S J . (originally 
following unified form: 



(S) 

tions ij.j- . (originally defined in ref. [1]), which appear in eq. (4.18) can be written in the 



P hlft = ( C f*t + C frt ~ C/ir)T7T- + ( C fir + C fit ~ C 'frt 



+ ( C fir + C frt ~ C fit) 



Sit^rt Sit Zt,ir 

1 Zr.it 



(E.2) 



Srt Zt ir 



This general form hides the fact that ft = g, which also implies Cf it = C/ 4 and Cf rt = Cf T 
(used in eq. (4.39)). Furthermore, according to our definition of the momentum fractions, 
we have 

Zi,jk UiQ Zij 



Zj,ik VjQ Zj^i 



(E.3) 
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for any i, j and k. Therefore, performing the integration over the soft phase space 
[dpi m+ i(pt] Q)] first, we find that the three terms in eq. (E.2) all lead to known func- 
tions. The eikonal term gives (—1) times the soft function J\(Y, e, yo, d' ), while both soft 
collinear terms give \ times the same soft-collinear function, /Ci(e, yo, d' ). For these func- 
tions we have the integral representations discussed in appendix C. Hence, we can express 
the integrated soft-triple collinear counterterms as linear combinations of two types of basic 
integrals, 

4 1} = ^T® 2e [ [dpQ\p r ,P i r-,Q)}^4jf(^,a,d(m,e))^]C 1 (e,y ,d' ), (E.4) 



and 



4 2 ) = _i|LQ2e /'[ d pg) ( ^^.g )] J_^ T/(a0)a)d(m)e))i7l(y? ;e , yo ,4). 

(E.5) 

The soft-collinear function 1C does not depend on any kinematic variable, therefore, 
it f actor izes completely from the integral over [dp^^]. The remaining integral can be 
expressed with the collinear function X(x\ e, cto, do — 1 + e, 0, &,0, 1) (the parameter do is 
shifted because this integral corresponds to m partons in the final state), so 

l£\x;e,a Q ,d Q ,y ,d ;k) = ^JCi(e,y ,d' )li(x;e,a ,d - l + e;k). (E.6) 

The soft integral depends on q, that has to be expressed with the integration 
variables of [dp^ ]: 

a i a + (1 — a ) x ^\ 

Y lr^ X Tr^^) = [a + {l _ a) ][a + {l _ a) {l _ v)] > 



where we used the usual abbreviations a = and v = v^~. Clearly, ^(x— ; a, v), as 

i r ir J ' i r ,Q\ ir ' ' ' ' 

well as the phase space measure in eq. (C.2) is symmetric in v <->■ 1 — v. Thus, it is sufficient 
to consider only the following coupled collinear and soft integral, 



l { g'(x;e,a ,do,yo,d' ;k) = -x / da (1 - a) 2a!o " 3+2e a" 1 " e [a + (1 - a)x] _1_£ 

Jo 

Jo \2a + (l-a)xJ 

x Ji(Y^ ? Q (x,a,v);e,y Q ,d'o) , k = -1,0, 1,2. 

In terms of the integrals ijp and Tg the integrated soft-triple collinear subtraction term 
can be written as in eq. (4.19). 
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Integrated soft-soft-collinear counterterm. To obtain the counterterm in eq. (4.20), 
we again perform the integration over the factorized soft phase space [dpf ' m , i(pt\ Q)] first. 
Now, we must distinguish two cases: (i) j, I and ir are three distinct labels, (ii) j or / coin- 
cide with ir. (We always have j 7^ I.) In the first case, the integral over the factorized soft 
phase space simply leads to the J\ function of eq. (C.10), and we find that the integrated 
soft-soft-collinear counterterm can be expressed as a linear combination of the integrals 

x (3) = i^! g2 e r [4p oj)] J_ z k /(ao , a _ )d(m , e)) 

^ J s ^f ' (E.9) 

x J 1 {Yj 7Q ;e,y ,d' ,0) . 
Using the collinear mapping formula p^ = (1 — )Pji> we fi n d 

Y-?n = ^^ = 7T \J 3 =Y-Jn- (E-10) 

b jQ b lQ V 1 a ir) b jQ b lQ 

Consequently, the collinear and soft integrals decouple, and we obtain 

X S ( x , Y 'i e > a o, d ,yo, d' ; k) =l 1 (x;e,a ,d - 1 + e; k) Ji(Y; e, y , d' ) . (E.ll) 

In case (ii) if e.g. j = (ir), the eikonal factor Sr ir )i(t) evaluates as in eq. (2.22), and 
the integral over the soft phase space leads to the one-mass soft function J^ lm \ defined in 
eq. (C.22) and computed in eq. (C.25). Hence we find that the integrated soft-soft-collinear 
counterterm can be expressed as a linear combination of the integrals 

4 4) = '^fQ 2 ' [ [dpg?]— 4 ? /(«o,« ?f ,d(m,6)) 

De J 8 ir ' (E.12) 

r(lm) 



xJ {lm) (Y Ci?)7tQ ,P { - ?) ;e,yo,d' ). 



This time the collinear integral does not decouple because the parameters Y^-p-u q an d 
P(7f\ depend on the collinear integration variables. Using the definition of pf r and pf, 
we have p'f., = (1 - a^ ? )pf r + a^-Q* 1 and = (1 - a^ ? )p(\ Consequently, 



i r \ i r liJ vr Q 

and 

j (x,y;e,a ,d ,?/o,do;A;) = x / da (1 - a) 2do_3+2e a" 1-e [a; + (1 - a)x] 

Jo 



(E.13) 



-1-e 



1 dvv-*(l-v)-<( a + ( 1 - a) ™) k (E.14) 
x ^ {lm) (^ ?) T iQ (^, y, a), % ?) (x, a); e, 2/0, do) , 

(A; = -1,0,1,2). 

In terms of the integrals and xj. , the integrated soft-soft-collinear counterterm 
can be expressed as in eq. (4.21). 
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Integrated soft-triple collinear-soft-collinear counterterm. When computing the 
integrals in eq. (4.22), we first perform the integral over the soft phase space. Using the 
same frame as in eq. (C.23), we find that the integral leads to the one-mass soft-collinear 
function JC^ lm ^ defined in eq. (C.33) and computed in eq. (C.36). Thus the integrated soft- 
triple collinear-soft-collinear counterterm can be expressed as linear combination of the 
integrals 

x (5) = 16^ g2e f [dpg)] J_4 /(ao>tt?p , dim ,e)) 

De J S ^r ' (E.15) 

The soft-collinear integral K,^ lm \f3) does not decouple from the collinear one because j3 = 
P(ir)^yirQ' a lr ) depends on the collinear integration variable as in eq. (E.13). Thus, 

, , rao 
l s '(x; e,a ,do,y , d' ; k) = x / da (1 - a) 2d °-' i+2e a- l ~ e [a + (1 - a)x] _1_e 

Jo 

x [ 1 dvv-<(l-v)-<( a + ()- a ^) k (E.16) 
Jo K ' \2a + {l-a)x) 

x /C( lm )(/3 (?f } (x, a); e, y ,d' ) , k = -1,0,1,2. 

(5) 

In terms of the integrals Tg , we find the result as in eq. (4.23). 
E.2 The soft-double soft-type counterterms 

The remaining soft-type integrated counterterms involve two successive soft momentum 
mappings, which leads to an exact factorisation of the original m + 2-particle phase space 
in the form 

d<£ m+2 (M) = d^ m ({p}( n ? '*))[dpg(p ? ;Q)][d^ +1 (p t ;Q)] . (E.17) 

The factorized phase space measures [dp[ r ^] and [dp^ m+1 ] are given in eq. (C.12), with 
appropriate changes in labelling (including m — > m — 1 in the first case) . We often need to 
express the two-particle scaled invariants of the momenta after the first mapping ('hatted 
momenta') with the final momenta obtained after the second mapping ('tilded momenta'). 
The relevant formulae, collected here for later reference, are: 

vii = i 1 - vrQ)y~u > yir=y~kri 

(E.18) 

ykQ = ( 1 -yrQ)y~ kQ + y~kr > k = i,j,i. 

Integrated soft-triple collinear-double soft counterterms. In computing the inte- 
grals in eq. (4.24), we first perform the integration over the soft phase space factor [dp^ , -J. 
This integration leads to either the J\ soft function, or the K,\ soft-collinear function of 



- 75 - 



section C. Then the remaining integral over [dp 1 ^] involves 

£ J if r A 



(E.19) 

x |/Ci(e,y )rfo), -^(^Tf^je^o^o)} ■ 

When the result of the first integration is a JC\ function, it decouples from the second 
integral, which gives a soft-collinear function again. Thus we find that all terms in the 
integrand of the type 

1 2 



Sjl Zljj S^J z j ti 

integrate to the product of two soft-collinear integrals: 



(E.20) 



^4 6) (e,yo,^o) = l;fci(e,yo,d' - 1 + e)/Ci(e, y Q , d' ) . (E.21) 

On the other hand, when the result of the first integration is a soft function J\, the 
two integrals do not decouple. In order to compute the second integral over the phase 
space factor [dp^], we choose the usual frame (C.ll), with orientation specified by 

= £ ? (1,...,1). (E.22) 

Using eq. (E.18), we compute 

Y ^ = Vj 9 = yj? = i-cos^ 

i9 ' Q VlQVrQ [{ l -yrQ)y~iQ+y~ i9 ]yvQ 2-y ?Q (l + COS?9) ' 



Then we find that the term in the integrand of the type 



(E.24) 



leads to the integral 

lP(e,y ,d> ) = -2 / dyy-^l-yfo-^ 
Jo 

x f dzz'^il - z)~ £ (l - y + yz) (E.25) 
Jo 

x Ji ( — ; e, y Q , d' ) . 

\l-y + yz J 

To write eq. (E.25) in the above form, we made the usual substitution of cost? — > 1 — 2z. 

In terms of the integrals Tg\e, yo, d' ) and Ig\e, yo, d' ) the integrated counterterm in 
eq. (4.24) can be expressed as in eq. (4.25). 
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Integrated soft-soft-collinear-double soft counterterm. In the case of the integral 
in eq. (4.26), we begin by integrating the eikonal function ^Sji(t) over the soft phase space 
[clpfL+J- However, whenever j or I is equal to ir (recall that j 7^ I), the eikonal factor 
evaluates as in eq. (2.22). Therefore, we have to distinguish two cases: (i) j, I and ir are 
three distinct labels, thus the integration over the first soft phase space leads to a J\ soft 
function and we obtain the integral 



r(8) 



167T 2 



Q 2e / [dp 



% r r , 1 



f{yo,yrQ,d'(m,e)) 



(E.26) 



4" = ^ JvrfSi 



1 Z'- 



it r , i 



and (ii) j or I coincide with (ir), hence the integration over the first soft phase space leads to 
a ,j( lm ) one-mass soft function and we find the integral (choosing j = ir for concreteness) 

.2 

-f(yo,y?Qid'(m,€)) 

(E.27) 

x^ (lm) ( y (?f)TQ'%.)^^o,d' ). 

To proceed, we must express the parameters of the soft functions with independent mo- 
menta. In the first case, using eq. (E.18), we can express Yyj g as 



Y 



ji,c 



4Y- 



ji,c 



[2(1 - y ?Q ) + 2y~. ? /y 7o ][2(l - y ?Q ) + 2y~ l? /y~ l 



(E.28) 



In the second case the one-mass soft function also depends on the velocity of the momentum 



(ir) 



(E.29) 



(ir)Q 



Again using eq. (E.18), we find 



2(1 - y?o) + 1y~ i? /y~ iQ + 2 y?Q /y~ iQ - I6y~ i? ly~ 



iQ 



(ir) 



Y, 



2(1 - y?Q) + 2y~ i9 ly~ iQ + 2y ?Q /y~ iQ 

4(1 - y? Q )Y TlQ + 4y l9 /(y~ iQ y lQ ) 



(ir)l,Q 



2(1 - y? q) + 2y~~jy~ iQ + 2y ?Q /y~ iQ 2(1 - y 9 Q ) + 2y J? /y~ u 



(E.30) 
(E.31) 



Turning to the integral over [dp^], we use two different orientations of frames in the 
two cases. In the first case, we fix the orientation such that 



P. 



5 " " " J 7 5 



1) 



P 



,smxy,cosx 



(E.32) 



that implies 
2y~.~ 



2/?q(1 — sinxj sin # cosy? — cos xj cost?) , (E.33) 
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where 

cos xj = cos x( Y jJ :Q ) (EM) 

with eq. (D.52) for cos%(y). Thus we see that although the integral over [dp^'^] is of 
soft-collinear-type, we cannot trivially integrate over ip because Yj^ q depends on it. Thus 

x j d(cosi?)d(cos^)(sin^)- 2e (sin^)- 1 - 2e (E.35) 

where the explicit dependence of the argument of the soft function (Y-rf n in eq. (E.28)) 
on the integration variables is 



-i 



(E.36) 



Y n>Q (Y, y, 0, ip) = AY (I - y) (2 - y(l + cos 0)) 

x ^2 — y (1 + sin x(X) sm $ cos V 3 + cos x(F) cos ^) 
In the second case, the orientation of the frame is fixed by 

Pi = E l( l i--- '!)' Pi = E l( 1 i ■■•>sinxj,cosxj), (E.37) 

that implies 

— — = 2/f o(l ~~ cos j — = 2/fo(l ~~ sm Xr sin cos (/3 — cos^t cost?) , (E.38) 

yiQ y~i Q 

where 

cos xj = cos x( Y jJ tQ ) (E.39) 
with eq. (D.52) for cosx(y). Then 

if (x, Y, e, y , d' ) = T dy y~^(l - 

x Z 1 d(cos 0) d(cos y) (sin fl)- 2e (sin v?)' 1 ^ 2 ~ ^ + c ° s ^) (E.40) 
,y_ 1 1 — cos v 

x J^(Y G9)lQ (x,Y,y^^),f] Ci?) (y^^);e,yo,d' ) , 

where and Y^^j q are given in eqs. (E.30) and (E.31), with the ratios of invariants 

in eq. (E.38) (x corresponds to y~Q, Y to Yrj q and y to j/f q). 

Our final result for the integrated counterterm defined in eq. (4.26) is presented in 
eq. (4.27). 
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Integrated soft-triple collinear-soft-collinear-double soft counterterm. The in- 
tegral over the soft phase space factor [dpf ' m , -J in eq. (4.28) is a one-mass soft-collinear 
integral JC^ lm ^ of eq. (C.33), which is computed in eq. (C.36), where the velocity of the 
momentum pf + p% is found in eq. (E.29). Then the integral over the measure [dp^], 

x (io) = i_^ Q2e r ^ _L^. f(yo>ytQ ^ M) 

S i r Z r,i (E.41) 



x ^ (lm) (%f); e ^o,^o), 



is again of soft-collinear-type, but the two integrals are coupled through that de- 

pends on the integration variables of the second integral as in eq. (E.30). The ratio of the 
invariants 2y~ i? /y~ iQ is 

= y ?Q (l-costf), (E.42) 

y iQ 

in a frame whose orientation is fixed by setting 

pf = E 1 (1,...,1). (E.43) 
Thus, [S^C^CS^S^] is equal to the integral 



4°\x;e,y ,d' ) = 2 dy y^(l - 



x / dzz-^il- z)~ e (l-y + yz) ( E - 44 ) 



x^ lm ^(3(x,y,z);e,y ,d' 
with x corresponding to y~ n , y to y^o and 



i -y + yz + y/x 

where, as usual, we set cosi? — > 1 — 2z. 

Our final result for the integrated counterterm defined in eq. (4.28) is presented in 
eq. (4.29). 

Integrated soft-double soft counterterms. There are two types of integrated soft- 
double soft counterterms: an 'abelian' one in eq. (4.30) and a 'non-abelian' one in eq. (4.31). 
Let us first consider the 'abelian' case. Performing the integral over [dpjj*^ , -J first, we 

obtain a soft function — Ji^Y^I qj e i Vo, d' ). In order to compute the integral over [dp^], 

4%i = [ [d&W k (?)f(yo,yr Q ,d'(m,e)) 

^ J 4 (E.46) 
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we must first express Yj-j- q in terms of Pj and pf as in eq. (E.28) (j, I ^ r), which is seen 

to depend on the integration variables in [dp^^] through the appearance of and yj~ 

in the denominator. Thus, the integration over [dp^J, which is of soft-type, is nontrivial. 
Also, although i ^ k and j 7^ I, there is no restriction on whether or not i,k is equal to 
j, I. Thus we must consider the following three cases: (i) all of i, k, j and I are distinct, 
(ii) only three of the four indices are distinct and, e.g. I = k, and (iii) only two indices are 
distinct and e.g. j = i and / = k. 

Case (i) requires at least four hard partons in the final state. Hence the corresponding 
integrated counterterm, [S t S^ t )]( l > fe )C?>0 with all labels distinct, does not enter a computation 
of two- or three-jet quantities, and we will not consider it in this paper. 

In case (ii), we have Yj-j g — > Yj^ g, and this is expressed with the independent 

momenta as in eq. (E.28), after a I — > k replacement. To evaluate the integral over [dp 1 ^] 
in eq. (E.46), we use the partial fraction identity (D.45) for the eikonal factor S^-j: ( r ) (with 
the substitutions j — > i and I — > k). Further, we restrict our attention to the case when 
there are precisely three hard partons in the final state. As discussed around eq. (D.54), 
this leads to a constrained kinematics for the three momenta pf, p£ and pj, and we take 
this into account below. It is convenient to introduce two different orientations of the frame 
eq. (C.ll) and integrate the two terms of the partial fraction in these different frames. In 
the first, we set 



Pi = E lO-i • • • J 1 ) ' Pk= E %( 1 '-- ■ ' sin Xik, cos Xik) , 
Pj = E j( 1 ' - ■ ■ '- sin Xij,c°sXij) , 



(E.47) 



so 



-^=y fQ (l-costf), (E.48) 
y iQ 

k r 



y j r 



UrQ(l - sin Xik sin 1? cos (p - cos Xik cos i?) , (E.49) 
VrCji)- + s ^ n Xij sin "i? cos (p — cosxij cosi?) , (E.50) 



3 Q 



where 



cosxife = cosxO^ q) > cosxij = cos xQ^q) , (E.51) 

with cosxQO given by eq. (D.52). In the second frame we exchange i and k, whose only 
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effect on the integrand is to interchange Xij an d Xjk- Thus we find 

x J d(cos??)d(cos^)(sin^)- 2e (sin^)- 1 - 2e (l-cos^)" 1 

x ^2 — cos-i? — sinx(Yi) sin $ cos <p — cos%(Yi) cosi?^ 
X «^( y 'i*,o( lr i» 1 2 J I/ J t? ) ¥');e,ift),dd > 
+ J (Vy £ q (Yi , Y 3 , y , i?, if ) ; e, y , d' 



(E.52) 



where 

Y 3%p (Yi,Y 2 ,y,0,<p) = 4Y 3 (l-y) 



x ^2 — y[l + sinx(li) sin # cos tp + cosx(Yi) cos0]^ (E.53) 

x (2 — y[l — sin x(Yz) sin •& cosip + cos x(Y2) cos $]\ 

Above, Y\ = Yr^ g, Yi = Y-r~-g and Y3 = Y~-^ g. The two terms in the square bracket 
in eq. (E.52) correspond to the two terms of the partial fraction in eq. (D.45). Their in- 
tegral representations are formally identical, only the kinematic invariants Y2 and Y3 are 
interchanged. 

In terms of the integral 2^ - k the integrated subtraction term defined in eq. (4.30) 
can be expressed as in eq. (4.32). 

Finally, case (iii) is obtained trivially from the previous one. Setting j = i implies 

Y-r-r n -> Y^-r n = , Y~. r „ -> Y~ r n , (E.54) 
therefore, depends only on Yr-^ g. Then integral (E.52) simplifies to 

(u) fv......j\- av. r^ 1 " 6 ) r° j., .,-i-2e/i 



x / d(cos^)d(cos^)(sin??)- 2e (sin^)- 1 ~ 2e (l-cos^)" 1 
■/-i (E.55) 

x ^2 — cos — sinx(Yi) sin i? cos ip — cos x(^i) cos $1 



and it replaces Ig^ljf, in eq. (4.32). Above 

F T £ (Fx, y, 0, if) =4Y 1 (l-y)[2- y[l + cos 0] " 



_ i (E.56) 
x ^2 — y[l + smx(Y) sin0cos <^ + cos x(Yi) cos 0]^ 
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Turning now to the 'non-abelian' contribution, eq. (4.31), we see that we need to con- 
sider a single new integral, 



r (12) _ 



Q 2e Jld P ^S~ % (?)f(yo,yrQ,d'(m,e)) 



(E.57) 



(obtained by performing the integral over [dp^ +1 ] as before). We must express Yj-~ q in 
terms of (see eq. (E.18)), 



Y- 



2 y~ir/yi Q 



rr ' Q y-ioVrQ V?Q2Q--V?Q) + 2Vi 9 /yii 



(E.58) 



We use the partial fraction identity (D.45) to write the eikonal factor as a sum of two 
terms, and choose two different orientations of the frame (C.ll) for each term. In the first 
one we set 

pi = £ 7 (1,...,1), p£ = E~ k (l,...,smx~ k ,cos X ~ k ), (E.59) 

thus 

2 V~i 



y~i Q 

where 



VrCji 1 ~ cos $) 



= UfQ^X — sin x k sin t? cos ip — cos x k cos i?) , (E.60) 



y 



kQ 



cosXfc =cos X {Y~ i ~ kQ ), 



(E.61) 



with eq. (D.52) for cos x(Y)- The second frame is obtained by the interchange i •<-)■ k, which 

(12) 

again implies change only in the argument of the soft function. Then the integral I$. ik 
equals 



-1+e 



x y d (cos tf)d (cosy?) (sin^)- 2£ (sin<^)- 1 - 2e (l - cosi?)" 1 



x ( 2 — cos $ — sin x(X) sm ^ cos V ~ cos xOO cos $ 
1 — cos 7? 



(E.62) 



Ji 



+ Ji 



2 — y(l + cos w) / 

(1 — sinx(^) sin # cos </? — cos xOO cos $ 
2 — y(l + sinx(^) sin t? cos 99 + cosx(^") cos??) 



;e,2/o,d' 



In terms of the integrals 2^^ ^ and Ig. ik , the integrated subtraction term defined in 
eq. (4.31) can be expressed as in eq. (4.33). 
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F Integrating the soft-collinear-type terms 

In this appendix, we discuss the integration of the soft-collinear-type counterterms of sec- 

(i) 

tion 4.3. In particular, we define and give an explicit integral representation of all 
functions (i = 1, . . . , 3). 

F.l Soft-collinear-triple collinear-type counterterms 

The soft-collinear-triple collinear-type counterterms involve a soft momentum mapping 
followed by a collinear one, which leads to an exact factorisation of the original m+2-particle 
phase space in the form of eq. (E.l). To evaluate the integrals of eqs. (4.34) and (4.35) 
over the factorized one-particle phase space measures in eq. (E.l), we first observe that the 
integral over the soft measure is a JC\ soft-collinear function, and the remaining the integral 
over the collinear measure can be expressed as linear combination of the integrals 

4? = ^Q 2e |[d^ ) ]^4 ? /(ao,a ?f ,rf(m,£))X: 1 (e,yo,d( ) ), (F.l) 

which is just twice the integral in eq. (E.4), computed in eq. (E.6), thus 

l^\x;e,a ,d ,y ,d' Q ;k) = 2l£\x;e,a ,d ,y ,d' ;k) . (F.2) 
F.2 Soft-collinear-double soft-type counterterms 

The soft-collinear-double soft-type counterterms involve two successive soft momentum 
mappings, which leads to an exact factorisation of the original m + 2-particle phase space 
in the form of eq. (E.17). To evaluate the integrals of eqs. (4.36), (4.37) and (4.38) over the 
factorized one-particle phase space measures in eq. (E.17), we first observe that the integral 
over the soft measure [dp^ +1 ] is again a /Ci soft-collinear function, and the remaining 
integral over the second soft measure contains either an eikonal factor or its collinear limit, 




i r r , i 



(F.3) 



(3) 

where in we recognise (twice) the function already defined in eq. (E.19), hence 

lg\e,y Q ,d! Q ) = 2lf{e,y Q ,d' Q ). (FA) 

(2) 

The second integral decouples from the first one also for 1^ because 1C\ is independent 
of the kinematics. Then, the final integral over [dp 1 ^] gives a soft function J\, so in this 



case, 



1§(Y; e, yo, d' ) = MY; e, y , d' - 1 + e)/Ci(e, y , d' ) , (F.5) 



where Y corresponds to Y~j n . 

3 f >W 



In terms of the functions 2^, (i = 1, 2, 3) we find the results given in eq. (4.39). 
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